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XTi ABSTRACT 

The observed radial and vertical metallicity distribution of old stars in the Milky Way disk provides 
a powerful constraint on the chemical enrichment and dynamical history of the disk. We present 

T-H the radial metallicity gradient, A[Fe/H]/Ai?, as a function of height above the plane, \Z\, using 

I 7010 main sequence turnoff stars observed by the Sloan Extension for Galactic Understanding and 

ILJ Exploration (SEGUE) survey. The sample consists of mostly old thin and thick disk stars, with a 

.J^ minimal contribution from the stellar halo, in the region 6 < i? < 16 kpc, 0.15 < \Z\ < 1.5 kpc. 

The data reveal that the radial metallicity gradient becomes flat at heights \Z\ > 1 kpc. The median 
metallicity at large \Z\ is consistent with the metallicities seen in outer disk open clusters, which 
exhibit a flat radial gradient at [Fe/H] ^ —0.5. We note that the outer disk clusters are also located 
at large \Z\; because the flat gradient extends to small R for our sample, there is some ambiguity 
in whether the observed trends for clusters are due to a change in i? or \Z\. We therefore stress the 
importance of considering both the radial and vertical directions when measuring spatial abundance 
trends in the disk. The flattening of the gradient at high \Z\ also has implications on thick disk 
formation scenarios, which predict different metallicity patterns in the thick disk. A flat gradient, 
such as we observe, is predicted by a turbulent disk at high redshift, but may also be consistent with 
radial migration, as long as mixing is strong. We test our analysis methods using a mock catalog 
based on the model of Schonrich & Binney, and we estimate our distance errors to be ^ 25%. We also 
show that we can properly correct for selection biases by assigning weights to our targets. 
Subject headings: Galaxy: abundances. Galaxy: disk, Galaxy: evolution, Galaxy: formation 
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1. INTRODUCTION 

1.1. Metallicity Gradients 

The spatial variation in the metallicity distribution of 
old stars in the Milky Way disk is linked to the formation 
and evolution of the Galaxy. The metallicity of stars at a 
particular place in the disk depends on the gas accretion 
rate, star formation history, and subsequent evolution at 
that location. For example, in the simplest picture of 
"inside-out" disk formation, low angular momentum gas 
falls to the center of the halo first, forming stars ear- 
lier and becoming c hem ically enrich ed much faster than 
the outer disk (e.g., [Larson.,, 1976; Matteucci & Francois| 



1989| [Chiappini et al.||1997t |Pra ntzos & Boissier 20U0). 
In this scenario, heavy element abundances decrease as 
a function of Galactocentric radius, i.e., the metallicity 
gradient is negativ e. The presen ce of radia l flows (e.g. 



Lacey fc F all 1985; |Goetz k Koe ppen 19921 |Portinanfc| 
Chiosi 2000; SpitoniT^ Matteucci 2011 ) and the nature of 



the early infailing gas (e.g., Cresci et al. (20To| )), however, 
have significant impacts on the chemical evolution of the 
disk and can generate a gradient that is weaker or even 
reversed compared to the simplest picture. Thus, obser- 
vations of the radial metallicity gradient of the disk are 
crucial to constraining chemical evolution models (e.g., 



Chiappini et al.|2001[[Cescutti et al.|2007||Magrim et al.| 



200^ 



The radial metallicity gradient of the Milky Way disk 
has been measured using a number of different tracers. 



between ~ —0.01 and —0.09 dex kpc~^ (e.g., Caputo 


et al.|2001||Friel et al.|2002 


Chen et al.l2003| Luck et al. 


2006| 


Lemasle et al.||2008| 


Sestito et al.||2008| IPediceUi 


et al. 


2009 


Luck & Lambert 


201 ip. These tracers rep- 



resent the composition ot the mterstellar gas at the time 
that they were formed, thus the wide variety of tracers 
studies probe the metallicity gradient at different times. 
This simple picture, however, can be complicated by pro- 
cesses that change the orbits of stars, such as dynamical 
heating from perturbations like spiral structure, molecu- 
lar clouds, or minor mergers, which can make gradients 
shallower or wash them out completely. 

In addition to the steepness of the gradient, non-linear 
features in the metallicity distribution provide further 
observational constraints. Some authors, for example, 
have noted the presence of a discontinuity in the radial 
metallicity gradient at a Galactocentric radius of i? ~ 10 
kpc, beyond which the metallicity gradient becomes shal- 
lower or flat (i.e., slope close to or equal to zero). This 



(e.g., 


Twarog et al."1997^'Yong et al.|2005 


Carraro et al.l 


2007; 
et al. 


Sestito et al. 2008 ) and < 
2002b(|Yong et ai.||2006 


Jepheids (e.g., j 
Pedicelli et al" 


\ndrievsky 
2009 j). 


Possible explanations tor this discontinuity include dy- 



namical int eractions within the disk (e.g., Andrievsky 



et al. 2004 ) and merger or accretion events in the outer 
disk (e.g., Yong et aL||2005 ). In the former scenario, the 
presence of the Galactic bar and spiral arms may influ- 
ence the star formation rate and flow of gas throughout 
the disk, affecting the amount of chemical enrichment 
that occurs at different radii. The latter scenario invokes 
star formation triggered by small accretion events, which 
is thought to explain the chemical abundances observed 
for outer disk open clusters (in particular the enhanced 



abundances of a and r-process elements), which differ 
from the general Galactic population. 

In this paper, we present the radial metallicity gradi- 
ent A[Fe/H]/Ai? of the Milky Way disk, as a function of 
height above the plane \Z\, using a sample of 7010 field 
stars from the Sloan Extension for Galacti c Understand- 
ing and Exploration (SEGUE; Yanny et al. 2009), part 



of the Sloan Digital Sky Survey (SPSS; [York et al.,2000^ . 
The sample covers Galactic coordinates b < R < w kpc, 
0.15 < \Z\ < 1.5 kpc, where R is the cylindrical Galac- 
tocentric radius and \Z\ is absolute distance from the 
plane. Our field stars are older than the open clusters 
and Cepheids used in previous gradient measurements in 
the literature and serve to extend the observations of the 
metallicity distribution of the disk to older tracers, which 
can provide constraints on the strength of the gradient 
at early times and how much radial mixing occurred. 

Additionally, we examine whether a discontinuity in 
the radial metallicity gradient exists in the old disk stars. 
While the distances derived for individual field stars are 
less accurate than for other tracers, our sample is suffi- 
ciently large to divide into bins of \Z\, which allows for 
an examination of the metallicity distribution in the disk 
as a function of both R and \Z\. The distinction between 
R and |^| is important, as many of the outer disk tracers 
in the literature are also located far from the midplane, 
and the question of whether the reported trends are a 
function of R, \Z\, or both, needs to be assessed. 

1.2. Thick Disk Formation 

This work is further motivated by the idea that the 
metallicity gradient of the old disk may be used as an 
observational constraint to distinguish between possible 
formation mechanisms for the thick disk. Traditionally, 
the Galactic disk can be thought of as the sum of two 
components: a thin disk of young metal-rich st a.rs and a 



thick d isk of ol der, more metal-poor stars (e.g., [G ilmore 



etal.|1 995; Chiba & Beers|2000||Bensby et al.|20(Mr 



et al.|2008| ). The ex istence of the Milky Way 'slEick disk 



vezic 



was fi rst noted by |Yoshii| ( 1982 ) and Gilmore & Reid 



([1983'), and thick disks with sirnilar kinematics, struc 



ture, and stellar populations have since been observed to 
be a common feature in nearby spiral galaxies (Dalca^ 
ton fc Bernstein|2002|[Yoachim fc Dalcanton|2005n200gr 



2008|). 



The ubiquity of thick disks in external galaxies, as well 
as the similarity of their properties, suggests that what- 
ever process is responsible for their existence is important 
in the formation and evolution of disk galaxies. Further- 
more, because the thick disk is old, the properties of its 
stars can be used as a "fossil record" of the disk's early 
formation. How stars end up in a thick disk, far from the 
plane of the Galactic disk, remains an open question. 

Several mechanisms for thick disk formation have been 
proposed, four of which are discussed below. Within the 
context of the hierarchical structure formation predicted 
by ACDM cosmology, a thick disk may arise through (1) 
the puffing up or vertical heating o f a pre-existing thin 



20081 


Read et al.| 


Puree 


1 et al. 2009 



2008; 



Kazantzidis ei al.[ [i>0O8, 200^ 
et al.„2011|); (2) the direct ac- 
lat merged with the 



Galaxy (Abadi et al. 2003); or (3) star formation in an 
early turbulent disk phase during a period of high gas 
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accre tion (e.g., Brook et al.||2004 2005 Bournaud et al 
2009 1. Even in the absence ot cosmologicai accretion, a 
thick disk may also arise through (4) rad ial migra tion of 
stellar orbits (e .g., Schonrich fc Biniiey] [2009a b; jLoeb 
man et al.|[201l| ). 



Each ot these scenarios is motivated by both theory 
and observations. Halo merger histories in cosmologi- 
cai N-body simulations suggest that the types of merg- 
ers required by scenarios 1 and 2 are common; |Stew- 



art et al. (2008) estimate that 70% of Milky Way-sized 



halos have experienced a 1:10 merger within the last 
10 Gyr. Streams in the halo of the Milky Way (e.g. 



Newberg et al.||2002| [Belokurov et al."2 007|) and other 

cz-Uelj 



galax ies (e.g., Ibata et~al.||2U01 "Iviartfncz-Uelgado et al 



2010 1 provide observational evidence of such accretion 
events. Other recent work, however, has emphasized 
the importance of smooth gas accretion in the growth 



of disk galaxies (Brooks et al. 2009 ; Dekel et al. |2009[ ), 



which can lead to a turbulent disk with high velocity dis 
persion as in scenario 3. Observational su pport for this 



picture include t he clump-cluster galaxies (Elmegreen & 
Elmeg rccn 2005) and the thick chain and spiral galax- 
ies (|Elmcgroen & El megrccn _2006 ) seen at hi gh red- 
shift^ Lastly, ISeTlwood &: Bmney ( 2002 ) and Roskar 



a) showed that resonant interactions be 



tween stars and transient spiral waves can change the 
radii of stellar orbits while keeping them on circular or- 
bits, leading to the kind of radial mixing necessary for 
scenario 4. Observations of nearby stars indicate that 
radial migration is an important process that may shape 
the correlations between the kinematics, metallicities, 



and ages of stars in the solar neighborhood (Haywood 



2008 

Wliile there is evidence that the mechanisms described 
above — minor mergers, early gas accretion, and radial 
migration — are at play in galaxy formation, the question 
of which mechanism, if any, is the dominant force behind 
thick disk formation remains unanswered. Recent work 



by Sales et al. ( 2009 ) is an example of how the kinematics 
predicted by the tour different scenarios — in particular, 
the distribution of stars' orbital eccentricities — can be 
used to test the various scen arios. This approach has 
been taken observationally by |Dierickx et al" (2010) and 
Wilson et al. (2011) using stars from t he SPSS and the 



RAdi al Velocity Experiment (RAVE; Steinmetz et al. 
20061, respectively. Both studies disfavor the minor 



merger scenarios, as they do not observe enough stars 
with high orbital eccentrici ties. Numerical simulations 
by ,Di Matteo et al. (|2011 ), however, suggest that the 



observed eccentricity distribution can be obtained in the 
minor merger scenario given different orbital parameters 



and satellite properties than those used in the Sales et al 
(|2009) analysis. The conflicting interpretations show the 
need for other observational constraints. In this paper 
we examine the radial metallicity gradient A[Fe/H]/Ai?, 
as a function of height above the Galactic plane \Z\, in 
order to further distinguish between different thick disk 
formation scenarios. 

For example, if thick disk stars originated in a thin disk 
and were subsequently heated by a minor merger as in 
scenario 1, then the observed radial metallicity gradient 
depends on the amount of mixing in the radial direction. 
If most of the heating occurs in the vertical direction, 
the thick disk will have the same metallicity gradient as 



the initial thin disk. Simulations, however, show that 
there can be substantial heating in the radial direction 



Hayashi & Chiba 2006 Kazantzidis et al. 2009 



Bird et al. 20lTj), which suggests that the gradient may 



be more shallow or flat than expected. If thick disk stars 
originated outside of the Galaxy and were deposited in 
the thick disk via accretion events as in scenario 2, then 
the metallicity distribution may ex hibit dumpiness. The 
simulations of 
disrupted sate 



Abadi et al. (2003) showed that a single 
lite roughly ends up in a torus of stars; 
several disrupted satellites would make up the thick disk 
by contributing stars of different metallicities at different 
radii. 

If thick disk stars originated in a turbulent gas disk 
at high redshift as in scenario 3, the short timescale for 
star formation makes the thick disk che mically homoge- 
neous, with no m etallicity gradient (Brook ct al. 2005j 
Bournaud et al. 2009). If thick disk stars originated in 
a thin disk and were pushed to larger radii through ra- 
dial migration as in scenario 4, then the original gradient 
would become washed out into a shallow or nonexistent 
radial metallici ty gra dient { Roskar et al.|2008a Sanchez- 
Blazquez et al.||2009 ). The presence ot a gradient, then, 
would rule out chemical homogeneity (scenario 3), while 
the strength of the gradient would constrain the amount 
of disk heating by minor mergers (scenario 1) and radial 
migration (scenario 4). Examining old disk stars in a 
large volume, beyond the solar neighborhood, will allow 
one to distinguish between the various scenarios. 

In contrast to many previous studies, we do not assign 
our stars to a thin disk or thick disk component. For 
samples of nearby stars in the solar neighborhood, this 
division is often done by assuming that the thick disk has 
a larger velocity dispersion and a slower mean rotation 
(i.e., th ick disk stars are k inem atically hot ) , as i n the 
cases of |Bensby et al.| ( |2003| and |Venn et al.| ( |2004| . An- 
other method ot se parating thin and thick disk stars is by 
their chemistry, as Lee et al. (2011a) do using [Fe/H] and 
[a/Fe]; they favor this type ot division because a star's 
composition is less li kely to change than its spat ial lo- 
cation or kinematics. Schonrich & Binney (2009b) have 
used mock observations to show that separating thin and 
thick disk stars using chemistry versus kinematics yields 
samples with different propertiesr] In addition, whether 
the thin and thick disks are truly distinct components 
is still an open question, with some studies arguing that 
the two components arise from a smooth correlation be- 
tw een chemical and ki nematic properties (see discussions 
b y|Haywood||2008| and [Ivezic et al.||2008| ) . 

Because our sample is not restricted to the solar neigh- 
borhood, we can compare the stellar populations of the 
thin and thick disks based on stars' locations instead of 
their kinematics or chemistry, which allows us to avoid 
assuming a specific model for the Milky Way disk. In 
this paper, we do not assign individual stars to the thin 
or thick disk. Instead, we will use the term thick disk 



^ For example, a star assigned to the thin disk using chemical 
criteria may be assigned to the thick disk using kinematic criteria 
if it belongs to the inner disk and is in the tail of the rotational 
velocity distribution. According to the| Schonrich fc Binney || |2009b|l 
model, this explains the tail of thick 3isk stars witn nigji [l^c/MJ 
and why most of the thick disk stars in the Bensby ct al. (20031 
and |Venn et al.| l |2004| samples are located at radii within the solar 
circfe 
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to refer to stars that are currently found at large dis- 
tances from the plane on their orbits; the term thin disk 
refers to stars that are found close to the Galactic mid- 
plane. Based on double-exponential fits to the vertical 
scale heights of the stellar de nsity distribu tion in the disk 
(e.g., [juric et al. 2008; dc Jong et al.|2010 ) we expect the 
thick disk to be the dominant population above \Z\ ^ 1.0 
kpc. 

Previous analyses of samples outside of the solar neigh- 
borhood have found no radia l metallicity gradient at ver 
tical heights \Z 



Latitude" pointings, which are restricted to Galactic lat 
itudes 8° < |6| < 16 



2009) 



> 1.0 kpc (Allende Prieto et al. 2006 



Juric et al.|2008). Our sample is complementary because 



our lines ot sight are located at relatively low Galactic 
latitude and we can directly compare the radial metallic- 
ity gradients of the thin and thick disks (up to \Z\ = 1.5 
kpc) using the same sample. In addition, we can explore 
whether the reported discontinuity in the outer disk is a 
purely radial trend or if a vertical trend is contributing 
to the observed flattening of the gradient at large R. 

The paper is organized as follows: The sample selec- 
tion and data are described in f|2] We then describe our 
methods of determining distances and correcting for the 
selection function in §S^and|4j respectively. Our gradi- 
ent measurements are presented in f|5j Error analysis is 
presented in ^ We discuss the results in S|7] and con- 
clude with a summary in fjsj For readers who are only 
interested in the results, we recommend skipping §S|4]and 
[6j Further description of our weighting scheme, intro- 
duced in f|4j is provided in the Appendix. Throughout 
our analysis we adopt the Galactocentric radius of the 
Sun, Rgc.0 — 8.0 kpc. 

2. DATA 
2.1. Sample Selection 
We measure the metallicity gradient of old main 



fields from the SEGUE survey (Yanm 


f et al. 


2009 


|Aihara et al.||2011| |Eisenstein et al.] 


2011 


Kockosi et al. 


2011, in preparation). These stars a 


How us to reach the 



largest distances probed by main sequence stars within 
the fixed magnitude limits of the survey . The data 
were obtained using the sarne tele scope (Gunn et al. 
2006|), camera ([Gunn et al.| 119981) , and filter system 



p-uFugita et al.|]l9y6| ) asTEe ^D5S. The old MSTO is 



selected to be m the blue part of the color-magnitude 
diagram (CMD), as described in detail below, and 
can be iden tified using the SDSS Data Release 7 
(DR7; |Abazaji_an et al. 20091 version of the Catalog 
Archive Server]^ as targets with sspParams . zbclass 
= STAR, SpecObjAll.primTarget = 2048, and 
PlateX.programName = seglow'/o. We also require 
that there are no repeat observations so that each 
star is only counted once. An equivalent qu ery 
in Data Release ^ (DR8; |Aihara et al.| |201l[ ) is 
SpecObj All. class = STAR, SpecDbjAll.primTarget 
= 2048, PlateX.programName = seglow"/,, and 
PlateX. isPrimary = 1, where the last requirement 
removes repeat observations. 

The programiName qualifier selects our targets from a 
subset of 22 SEGUE plug-plates that comprise the "Low- 

^ Ihttp: / /casjobs. sdss.org/CasJobs/ 

* http: / /skyservicc. pha.jhu.edu/casjobs/ 



see §3.15 of Yanny et al 

These lines of sight are high enough to avoid the young 
star-forming disk, as well as the regions with the most 
crowding and highest reddening, but also sufficiently low 
that they have a long sightline through the disk. The 
lines of sight fall into roughly two groups in Galactic 
longitude: seven at 50° < I < 110° and another four 
toward the anticenter, 170° < / < 210°. 

Each plate covers 7 square degrees on the sky with 
targets in the magnitude range 16 < g < 20, where the 
magnitudes have not been corrected for extinction. In 
this paper, we will refer to any reddening- or extinction- 
corrected magnitudes and colors with subscripts gsFP 
and go, f or correc tions derived from'Schlegel et al. (1998 



hereafter SFD98 ) and isochrone fitting (see jpl), respec- 
tively, lable [T] lists the properties of the II lines of 
sight (two plates per pointing) included in our sample, 
orde red by t he median extinction E{B — V), obtained 
from SFD98 For the total sample, E{B — V) varies be- 
tween 0.05 and 1.07 magnitudes. On average, there were 
600-700 spectra obtained per line of sight. 

For most of the SEGUE survey, which was at high 
Galactic latitude, targets are identified as MSTO stars 
based on their (u — g)sFi) and {g — r)sFD colors (see 



Yanny et al. 2009). In ugr color space, it is possible to 
separate the IvlS'l'O stars from metal-poor halo stars be- 
cause of the large ultraviolet excess of metal-poor stars. 
In the low latitude pointings, however, two issues arise. 
First, the w-band magnitudes and their uncertainties are 
unreliable due to the large extinction in these regions. 
Second, it is impossible to use a single constant selection 
in (g — r')sFD that will yield the same stellar population 
in every line of sight because the reddening in these fields 
is, on average, much higher and more variable than in the 
high latitude fields. For these low latitude fields, we use 
a targeting procedure that is more robust to reddening, 
which will reliably choose the stars at the blue edge of 
the CMD. Starting with the photometric objects identi- 
fied as stars in the imaging, this selection procedure is as 
follows: 

1. We remove all stars with g > 20 and i < 14.2 (using 
magnitudes uncorrected for extinction) to ensure 
that the targets are bright enough for high-quality 
spectroscopy in the expected exposure time. 

2. The 7 square degree area of each plate is large 
enough that the extinction can be highly variable 
across the plate, and there are always many more 
targets than fibers available. We remove the re- 
gions of highest extinction from consideration to 
maximize the number of useful spectra. For each 
half of the plate, we calculate the 75th percentile 
of the E{B — V) distribution using the total line 
of sight extinction from SFD98i This procedure is 
done for each half independently to ensure that the 
targets are approximately evenly distributed over 
the plates, since the reach of the fibers is only about 
half of the plate diameter. 

3. We remove all objects with E{B — V) larger than 
the higher of the two 75th percentile values. Tak- 
ing the higher value ensures that there are enough 
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TABLE 1 
Properties for 11 Lines of Sight. 



Plates RA(°) DEC {°) I (°) b (°) E(B - V)-^ N, 



2712 


2727 


105.6 


12.4 


203.0 


8.0 


0.09 


744 


2536 


2544 


286.7 


39.1 


70.0 


14.0 


0.15 


661 


2534 


2542 


277.6 


21.3 


50.0 


14.0 


0.17 


672 


2554 


2564 


303.0 


60.0 


94.0 


14.0 


0.19 


734 


2678 


2696 


98.1 


26.7 


187.0 


8.0 


0.24 


766 


2556 


2566 


330.2 


45.1 


94.0 


-8.0 


0.31 


728 


2668 


2672 


79.5 


16.6 


187.0 


-12.0 


0.33 


830 


2681 


2699 


71.5 


22.0 


178.0 


-15.0 


0.41 


758 


2537 


2545 


334.2 


69.4 


110.0 


10.5 


0.49 


708 


2538 


2546 


323.1 


73.6 


110.0 


16.0 


0.65 


716 


2555 


2565 


312.4 


56.6 


94.0 


8.0 


0.82 


511 



Median value for 



spectra in lino of sight using values from |SFD9"8l 



usable targets on each half of the plate to fill all 
the fibers given their limited reach across the plate. 
This should not bias the sample, as we do not ex- 
pect that the objects behind more extinction are 
intrinsically any different from those that are un- 
obscured. This is especially true for distant objects 
that are located far behind the dust. The magni- 
tudes used through out the r est of this procedure are 
correcte d using the \SFD9^ extinction values. The 
|SFD98| extinction was applied so that stars in the 
same approximate luminosity range were targeted 
in each line of sight despite the large variation in 
extinction among the different lines of sight. 

4. We examine the {g — r)sFD distribution in bins of 
5SFD, each one magnitude wide. For each distri- 
bution, we find the peak, which is the {g — r)sFD 
color of the MSTO in each ^-magnitude bin. In 
addition, we determine the half maximum on the 
blue side; this is [g - r)hait-max- 

5. The red cut for each bin is defined as {g — r)cut = 
{g - ?')haif-max + 0.25. We fit a fine to [g - r)cut as 
a function of (?bin, where gbin is the mean (7sfd of 
all the stars in each bin. 

6. All stars on the blue side of the line are defined as 
candidate spectroscopic targets with equal proba- 
bility of being selected. Targets are randomly cho- 
sen from the resulting candidate list. 

Though the (g — ?')sfd color of the population may 
change from field to field because of varying amounts of 
extinction, the identification of the MSTO stars as the 
bluer population holds in all lines of sight. As a result, 
this method is more robust to reddening than the stan- 
dard color cuts for normal SEGUE plates. Halo contam- 
ination is expected to be low in these plates, compared 
wit h th e higher latitude pointings (see further discussion 
in S6.2). The color cut, however, will bias our sample 
against metal-rich stars, which have redder colors. The 
severity of the bias depends on how much the MSTO 
color changes with metallicity, which in turn depends on 
the ages of s tars at eac h metallicity (see further discus- 
sion in §f|6T|and 6.3.2). 

Figure [T] shows the results of the procedure outlined 
above for two lines of sight with low and high extinc- 
tion (median E{B — V) = 0.17 and 0.41, respectively). 
The density of objects identified as stars in the photom- 
etry are plotted in grayscale and contours, while MSTO 



TABLE 2 

Cluster Metallicities Measured by the SSPP. 



Cluster'^ 


[Fo/H]Lit 


[Fc/H]sspp 


[Fe/H]Lit-[Fe/H]sspp 


M3 


-1.50 


-1.43 


-0.07 


M71 


-0.82 


-0.74 


-0.08 


NGC2158 


-0.25 


-0.29 


+0.04 


NGC2420 


-0.20'' 


-0.31 


+0.11 


M35 


-0.16 


-0.24 


+0.08 


M67 


+0.05== 


+0.00 


+0.05 


NGC6791 


+0.30 


+0.29 


+0.01 



^ Cluster data are also presented in |Smolinski et al.| ( |2bll[ | , with the 
exceptions of NGC2420 and M67. A full comparison of the SSPP for 
all data in the SEGUE cluster samples will be presented in Rockosi et 
al. (2012. in prep). 
Jacobson et al. KOllk. 

° psj!cTet ai.fproy . 



stars in our spectroscopic sample are plotted as blue cir- 
cles. All photometric objects bluer than the red limit 
of the spectroscopic sample were considered as candi- 
dates for spectroscopy, but only a randomly selected sub- 
set of those were actually observed. In the 22 low lat- 
itude plates, spectra of 7828 MSTO stars were taken. 
We keep the targets with good photometry and spectra 
with S/N > 10 per pixel, where each pixel corresponds 
to ^ lA. We remove any stars which have large dis- 
crepancies (> 0.1 mag) in f; — r color between different 
photometric reductions. We also remove rare blue stars 
with {g — r)sFD < —0.25. The resulting sample contains 
7655 spectra with a mean S/N ^ 30 per pixel. 



2.2. 



SEGUE Stellar Parameter Pipeline: Accuracy in 
Regions of High Extinction 



The SEGUE Stellar Parameter Pipeline (SSPP; Lee 



et al. 2008a) estimates the effective temperature Jeff, 
surface gravity log g, and metallicity [Fe/H] for each 
spectroscopic target in the survey. We use stellar pa- 
rameters from the version of the SSPP used for Data 
Release 8, which includes impro ved [Fe/H] estimates a t 
both high and low metallicities ( Smolinski et al.||2011 ). 
The SSPP has been extensively tested using globular and 
open clusters, where true cluster members are identified 
using t heir metallicities and radial velocities ( Lee et al. 
,2008b{ IS^linski et al.|[20Tl] ). 

We verity that these results hold for cluster members 
in our temperature and surface gravity range (5000 < 
Toff < 7000 K, log g > 3.3). Table [2j shows the com- 
parison between the literature values (column 2) and the 
SSPP for cluster members in our temperature and sur- 
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Fig. 1. — Main sequence turnoff (MSTO) selection. Spectroscopic targets (blue circles) are a randomly-selected subset of the MSTO 
stars identified using the photometry of all objects in the field (grayscale and contours). The contour labels indicate the number of stars 
per 0.375 by 0.1 magnitude q aFn-(q — ?')sFD bin. The results are shown for two lines of sight with low and high extinction (left and right, 
respectively), as measured by |SFD98| 



face gravity range (columns 3-4). The offsets in [Fe/H] 
between the hterature values and the SSPP are small 
(within 0.1 dex), and we see no trends with Tog, [Fe/H], 
and S/N. These tests show that we can reliably measure 
trends in [Fe/H] throughout our entire sample volume, 
and that the absolute values of the metallicities presented 
in the paper are accurate to 0.1 dex or better. 

Each parameter is estimated using multiple methods: 
11 for Teff, 10 for \ogg, and 12 for [Fe/H]. For T^e in 
particular, these methods include spectral fitting and x^- 
minimization using grids of synthetic spectra (kil3, k24, 
NGSl), measuring line indices (WBG, HA24, HD24), 
neural networks using training sets of both real and syn- 
thetic spectra (ANNRR, ANNSR), and q - r color pre- 



dictions (Tk, Tq, Ti). See Lee et al. (2008a) for complete 
details on each of these individual methods. Some meth- 
ods use only the spectra as input, some use only the pho- 
tometry, and some use both. The individual estimates of 
each parameter are averaged to obtain a final adopted 
value. Each individual estimate is valid for some range 
of (g — r)sFD color and S/N, which determines whether 
it will be included in the final average. A more compli- 
cated decision tree is used for [Fe/H] and is describ ed in 
detail in the Appendix of Smolinski et al. (2011). We 
have 7605 stars that have good stellar parameters; their 
temperatures fall in the range 5000 < Tcs < 7000 K, 
making them F and G dwarfs. 

Before using the SSPP parameters, we test whether 
the parameter estimates are affected by high Galactic 
extinction. The SSPP was designed to analyze the nor- 
mal SDSS and SEGUE data at high Galactic latitude 
and uses photometry that has be en correc ted for extinc- 
tion using the reddening maps of SFD98 These extinc- 
tion values reflect the total line-of-sight extinction, which 
means that the colors of less distant stars will be over- 
corrected; they will be too blue. This effect is likely to 
result in a systematic error in the parameters estimated 
by those methods in the SSPP that use the photometry. 

If this overcorrection affects our sample, we expect es- 
timates that use the photometry to be systematically 
different when the extinction is high. The left panels 
of Figure [2] show the discrepancy between the individ- 



ual Toff estimates and the adopted value as a function 
of E{B — V). Whether the estimate includes the pho- 
tometry is indicated in the bottom left of each righthand 
panel. We have only plotted estimates that were ac- 
cepted by the SSPP (i.e., the target falls in the (5 — r)sFD 
or S/N range in which the particular method is reliable). 
As expected, the photometry-dependent estimates of Toff 
(k24, WBG, Tk,Tg,Ti) are systematically higher (i.e., 
the color is bluer), compared to the adopted value, for 
the highest values of extinction. 

To remove this effect, we calculate new averaged values 
of Tcff using only the estimates from methods that ex- 
clude the photometry (kil3, ANNSR, ANNRR, NGSl, 
HA24, HD24). The right panels of Figure [5] show the 
discrepancy between the individual Tes estimates and 
the newly calculated spectra-only value as a function 
of E{B ~ V). In contrast to the left panels, there is 
no clear trend in AToff with extinction for the individ- 
ual spectra-only estimates, which suggests that our new 
spectra-only average of the temperature is more reliable. 
At E{B — V) > 0.8, the median discrepancy between 
the new spectra-only estimate and the original adopted 
value is about 100 — 200K, which is compara ble to the 
expected errors of the SSPP Tcff. We show in { 6.3.1 that 
this could amount to a systematic error in the distance 
of ~ 20 — 25%. For the remainder of the paper, T^s will 
refer to the spectra-only value. 

No trend with extinction is observed in the 
photometry-dependent estimates of [Fe/H] (k24, WBG, 
CaIIK2, CaIIK3, ACF, Call), so we keep the adopted 
values. This provides a more robust result, as the 
adopted value is an average of a larger number of es- 
timators. 

3. DISTANCES 

To calculate distances to each target, we use the 
spectra-only Teff and SSPP [Fe/H], plus the theoretical 
isochrones of An et al. (2009j^ which have been shown 
to be good matches to ugriz cluster fiducials. We do not 
make use of the SSPP log g estimate, as it is relatively 
inaccurate near the turnoff, where the expected range of 

^ http:/ /www. astronomy.ohio-statc.edu/iso/sdss. html 
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Fig. 2. — Left panels: Differences between individual SSPP temperature estimates and the adopted values as a function of 'SFD98| 
extinction E{B — V). For estimates that include photometry (k24, WBG, Tk,Tg,Ti), the deviation from the adopted value increases 
with extinction, indicating that the temperatures are overestimated due to overcorrection when using the SFD98 extinction values. Right 
panels: Differences between Individual SSPP temperature estimates and the spectra-only values as a function of extinction. The trend is 
no longer evident in the spectra-only estimates, indicating that the spectra-only temperature is more reliable for highly extincted objects. 
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gravities is smal l compared to the errors. Twelve sets of 
An et al. (2009) isochrones at metallicities in the range 
—3. (J < [i'e/HJ < +0.4 are available, with each set having 
a range of ages up to 15.8 Gyr. The distance uncertain- 
ties are discussed in §6.3. 1[ 

We assign all the stars in the sample to the isochrone 
with the closest metallicity. We then find the mean tem- 
perature for the stars in each metallicity bin by fitting 
a Gaussian to the distribution of effective temperatures. 
We identify the age of the isochrone with the turnoff 
temperature closest to the measured mean temperature; 
we refer to this as the turnoff age of the mean tempera- 
ture, or TAMT, for each metallicity. Targets hotter than 
the mean temperature cannot be placed on the TAMT 
isochrone, so younger isochrones must be used; these tar- 
gets are assigned to the oldest possible age (i.e., the old- 
est isochrone where the target is cooler than the turnoff). 

Figure [3] shows a schematic picture of how ages are 
assigned to targets depending on their effective tempera- 
tures. The top panel shows six representative isochrones 
at solar metallicity, while the bottom three panels show 
how stars are assigned to these isochrones based on their 
place in the Toff distribution. The procedure described 
above is shown in the panel labeled "TAMT" while the 
panels labeled "TO" and "ZAMS" show two other age 
a ssump tions that we use to estimate our distance errors 
(j |6.3.1| ). Table |3] lists the TAMT determined by finding 
the mean temperature for each metallicity bin for the 
7605 stars in our sample. 

Once we have assigned a given target to an isochrone 
with a particular age and metallicity, we use the 
isochrone to obtain the predicted g — r color for the tar- 
get's spectra-only Tcff by linearly interpolating on the 
isochrone in temperature-color space. A comparison of 
the predicted and observed g — r gives an estimate of the 
isochrone extinction in g — r, which is also used to deter- 
mine the extinction in the ff-band . The is ochrone extinc- 
tion is an improvement over the SFD98 values because 
it does not assume that the target lies behind all of the 
dust in the line of sight. We step through this procedure 
for the two isochrones with the nearest values of [Fe/H] 
and then linearly interpolate to find the predicted appar- 
ent and absolute magnitudes. The apparent magnitude, 
now corrected using the isochrone extinction, along with 
the predicted absolute magnitude in the g-band yield the 
distance. 

In addition to being used in the distance calculation, 
the isochrone extinction provides information about the 
dust distribution along different lines of sight in the field. 
The left panel of Figure [4] shows the isochrone extinction 
as a function of the derived distance for one line of sight, 
with the color indicating the SFD98 value. There is good 
general agreement between the two extinction estimates, 
especially on a relative scale. Quantitatively, the agree- 
ment is poorest for small distances where the targets are 
in front of some of the dust, but asymptotes to better 
agreement at l arge dis tances. This is consistent with the 
idea that the ISFD98I values are overestimates because 
they include all of the dust in the given line of sight. For 
the line of sight shown, most of the dust is located within 
2 kpc of the Sun, beyond which the isochrone extinction 
is approximately constant as a function of distance. 

The scatter in the left panel of Figure [4] may be ex- 
plained by the patchiness of the extinction on the plane 



[Fe/H] = 0.0 

1.12 Gyr 
1.59 Gyr 
2.51 Gyr 
5.01 Gyr 

8_9I Gvr 
15.9 Gyr 




Fig. 3. — Determination of isochrone ages for solar metallicity 
targets (—0.05 < [Fe/H] < 0.05, gray dots). Isochrones are shown 
in the top panel, with the temperature distribution in the bottom 
three panels. In practice, the set of isochrones used at each metal- 
licity includes all available ages, but for clarity we show only six — 
the oldest, youngest, and those located at the mean and 1-, 2-, and 
3-cr values of the temperature distribution. The colors in the lower 
panels indicate which isochrones in the top panel are used to cal- 
culate distances for those targets. We show how ages are assigned 
to our targets making three possible age assumptions: Turnoff Age 
of the Mean Temperature (TAMT), Turnoff (TO), and Zero-Age 
Main Sequence (ZAMS). TAMT: For each metallicity, most tar- 
gets are assigned to a single age (the TAMT) at which the turnoff 
temperature is closest to the measured mean of the temperature 
distribution of the SEGUE targets. Targets hotter than the mean 
are assigned the oldest possible age (i.e., the oldest isochrone where 
the target is cooler than the turnoff). TO: All targets are assigned 
to the oldest possible age, even those cooler than the mean of the 
temperature distribution. ZAMS: All targets are assigned to the 
youngest possible age. The TAMT assumption is used in our anal- 
ysis, while the TO and ZAMS assumptions are used to test the 
accuracy of our distance estimates (see i ]6.3.1[ l. 

of the sky. The right panel of Figure |4] shows that the 
distribution of E{B — V) varies on small scales. The 
scatter in the left panel is expected if each region of high 
extinction has a different dust distribution along the line 
of sight. The agreement between the isochrone extinc- 
tion and the |SFD98| values provides a sanity check which 
indicates that our isochrone extinction estimates are re- 
liable. 

Sixty-two stars in our sample (0.8%) end up with neg- 
ative values of isochrone extinction. Three kinds of stars 
fall into this category: (1) 17 stars are outliers with very 
blue g — r colors. These stars do not contribute to the 
final measurement; after we apply our weighting scheme 
they receive a CMD weight of zero (see Q . (2) 7 stars 
show large changes in g — r color between different pho- 
tometric reductions and are likely to be catastrophic er- 
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TABLE 3 

Mean Isochhone Ages as a Function of [Fe/H]. 



[Fc/H] 




[a/Fc] 


log (Totf) 


(Teff) (K) 


log (age) 


{age) (Gyr) 


Nobj 


[Fc/H] < 


-2.50 


+0.4 


3.768 


5864 


10.2 


15.8 


9 


-2.50 < [Fc/H] < 


-0.75 


+0.3 


3.768 


5864 


10.20 


15.8 


880 


-0.75 < [Fc/H] < 


-0.40 


+0.2 


3.775 


5958 


10.05 


11.2 


2695 


-0.40 < [Fc/H] < 


-0.25 


0.0 


3.772 


5912 


10.05 


11.2 


1527 


-0.25 < [Fc/H] < 


-0.15 


0.0 


3.772 


5912 


10.00 


10.0 


857 


-0.15 < [Fc/H] < 


-0.05 


0.0 


3.772 


5912 


9.95 


8.9 


732 


-0.05 < [Fc/H] < 


+0.05 


0.0 


3.763 


5796 


10.00 


10.0 


480 


+0.05 < [Fc/H] < 


+0.15 


0.0 


3.763 


5796 


9.95 


8.9 


229 


-1-0.15 < [Fc/H] < 


+0.30 


0.0 


3.763 


5796 


9.85 


7.1 


142 


-1-0.30 < [Fc/H] 




0.0 


3.763 


5796 


9.70 


5.0 


54 




Fig. 4. — Estimates of the extinction derived from isochrone fitting give a picture of the dust distribution along a given Une of sight. The 
left panel shows the isochrone extinction as a function of derived distance d from the Sun. At small d the isochrone extinction is smaller 
than the SFD98 value (indicated by symbol color), but at large d it asymptotes to the SFD98 value. For this line of sight, which has the 
highest median E(B — V) of our 11 fields, most of the dust is found within 2 kpc of the bun. The scatter in isochrone extinction for a 
given value of SFD98 extinction may reflect the patchiness of the dust on the plane of the sky. This effect is seen in the right panel, which 
shows multiple separate patches that have the same [SFD98| extinction values, but may have different distributions of dust along the line of 
sight. The color coding in the right panel is the same as in the left. 



rors, possibly due to blending in the relatively crowded, 
low-latitude fields. (3) The remaining 38 stars tend to 
be faint, and likely have negative reddening because the 
errors in their temperatures cause their predicted g — r 
colors to be redder than their observed colors. We only 
see these stars in fields where the reddening is low; in 
fields where reddening is higher, the observed g — r is 
much redder than the predicted g — r, so the tempera- 
ture errors are not large enough to cause stars to have 
negative reddening. We expect just as many objects to 
have temperature errors that cause their predicted g — r 
colors to be bluer, and we leave all of these objects in 
our sample, which should give us a more symmetric dis- 
tribution of errors in our distances. 

With distances for all of our targets, we obtain the spa- 
tial distribution of our sample in Galactic coordinates R 
and \Z\, as shown in Figure |5| The SEGUE targets are 
shown as blue dots. Our radial coverage extends mostly 
to Galactocentric radii at the solar circle and beyond, 
though at low \Z\ it is confined just to the solar cir- 
cle. About one-third of our sample is located below the 
Galactic plane; the 5 < 0° lines of sight cover the entire 
radial range outside of the solar circle. Thirty-three open 
clusters and 190 Cepheids analyzed in the literature us- 
ing high-resolution spectra are also shown (open symbols 



3.0 r — ' — ' — ' — ' — I — ' — ' — ' — ' — I — ' — ' — ' — ' — r 




5 10 15 20 

R (kpc) 



Fig. 5. — Spatial distribution of our sample (blue dots) in Galac- 
tic coordinates R and \Z\. We divide our sample into four \Z\ slices, 
indicated by the black boxes. Open cluster and Cepheid data are 
plott ed with open symbols and plusses/crosses, respectively (see 
j ]7.2| for details). Objects from the literature located at \Z\ > 1.5 
kpc are plotted in gray. We note that all of the outer disk clus- 
ters from the literature are located at large |Z|. Multiple literature 
values for a single cluster are connected by lines. 

and crosses, respectively). These data have been used to 
study the radial metallicity gradient of the disk and will 
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serve as comparison samples. The properties of the open 
clusters are listed in T able [4} and the data are described 



in more detail in [7.2 

We use our sample of MSTO stars to measure the 
metallicity gradient of the old disk in four different \Z\ 
slices as indicated by the black lines in Figure [5j Taking 
thin and thick disk scale heights to be ^ 300 and 900 
pc, respectively, the lower two slices (0.15 < \Z\ < 0.25 
and 0.25 < \Z\ < 0.5 kpc) are dominated by the thin 
disk. The third slice (0.5 < \Z\ < 1.0 kpc) is made up 
of a mix of the thin and thick disks, and the fourth slice 
(1.0 < \Z\ < 1.5 kpc) is dominated by the thick disk. Of 
our sample, 7180 stars fall into these four slices. 

Dividing our sample in this way allows a comparison 
of the radial metallicity gradient of the thin disk to that 
of the thick disk, as well as to distinguish between ra- 
dial and vertical trends. We note that all of the distant 
clusters that have been used to study the behavior of the 
radial metallicity gradient in the outer disk (i? > 15 kpc) 
are located at least 1.5 kpc from the Galactic midplane; 
these are shown in gray in Figure [5j 

4. CORRECTING FOR SELECTION BIASES: WEIGHTS 

To use field stars to determine the metallicity distri- 
bution in the disk, we must understand how the spectro- 
scopic sample is drawn from the underlying population. 
As there are many more stars than it was possible for 
SEGUE to obtain spectra, we must assess whether the 
spectroscopic sample is truly representative of all of the 
stars in the disk. It likely is not, because our stars are 
selected to be the bluer stars in the CMD, making our 
selection biased against metal-rich and older stars, which 
have redder colors. The severity of the metallicity bias 
will depend on the ages of metal-rich stars; the older and 
more metal-rich they are, the more likely they are to fall 
out of our sample. To correct for this bias, we employ 
a weighting scheme to step backward in our sample se- 
lection and reconstruct the properties of the underlying 
parent population. We describe the scheme briefly below; 
further details are in the Appendix. 

There are three major ways in which the spectroscopic 
sample is different from the parent population along each 
line of sight: (1) The photometric objects in regions 
with the highest extinction were not considered for spec- 
troscopy. (2) Not all candidates for spectroscopy are ob- 
served. (3) We observe only MSTO stars using a color 
cut that is biased against redder metal-rich stars. 

Each star in our sample is given three weights corre- 
sponding to the three differences listed above: (1) the 
area weight Wa, which depends on the coverage of tar- 
gets on the plane of the sky in each line of sight; (2) the 
CMD weight M^cmd, which depends on the target's loca- 
tion in the CMD and corrects for the random selection of 
a subsample of all candidates for spectroscopy that pass 
the MSTO selection; and (3) the LF weight Wlf, which 
depends on the target's Tes, [Fe/H], and location in the 
CMD, and corrects for the metallicity bias of the MSTO 
selection. The total weight W is the product of the three 
weights Wa, W^cmd, and Wlf- After removing targets 
with = 0, we are left with a sample of 7010 stars. De- 
tails about how each of these weights is calculated can 
be found in the Appendix. 

Figure |6l shows the distribution of the total weight 
(black) and each individual weight, as a function of R, 
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Fig. 6. — Variation in weights as a function of _R and \Z\. The 
variation in the total weight W (black) is mostly dependent on the 
CMD weight Wcmd (blue). The LF weight Wlf (red) is fairly 
constant at all locations, with the width of the distribution be- 
ing ~ 0.5 — 1.0 dex. The area weight Wa (green) is the smallest 
contribution and does not vary greatly between different lines of 
sight. While WcMD shows the most dramatic variation, it is also 
less uncertain than VFlf because it only requires counting objects 
in CMD bins. 

in four slices of \Z\. The panels on the right show the 
distribution of weights. Wa (green) is the smallest con- 
tribution and does not vary significantly. Wlf (red) is 
relatively constant as a function of R and \Z\, with the 
distribution having a width of ^ 0.5 — 1.0 dex. Since M^lf 
is relatively constant, a systematic error in VFlf, which 
could arise from using the wrong luminosity function, will 
not cause a spurious change in the ratio of metal-poor to 
metal-rich stars. 

WcMD (blue) shows the most variation because it nor- 
malizes for the fact that there are more stars in the in- 
ner disk than the outer disk. Although the variations in 
WcMD are large, this weight is less uncertain than Wlf^ 
as it only requires counting objects in each bin. Figure 
shows that the change in the total weight W mostly fo" 
lows the change in Wcmd- 

5. RESULTS: RADIAL METALLICITY GRADIENTS 

The total weights are applied to each target, allowing 
us to use our sample of MSTO stars to estimate the prop- 
erties of the underlying parent population. We divide our 
sample into four slices of \Z\: (1) 0.15 < \Z\ < 0.25, 
(2) 0.25 < \Z\ < 0.5, (3) 0.5 < \Z\ < 1.0, and (4) 
1.0 < \Z\ < 1.5 kpc. Within each \Z\ slice, we fit a 
linear gradient to the data, weighting each target by the 
total weight determined using the scheme described in 

Figure [71 shows the radial metallicity gradients for all 
four \Z\ slices for both the unweighted and weighted 
cases (blue triangles and red squares, respectively). The 
weighted slopes are written in the bottom left corner of 
each panel, with the unweighted values in parentheses. 
The quoted errors are derived from the Monte Carlo sim- 
ulations described in ^6.3.2 and include only the random 
errors from uncertainties in the stellar parameters. See 
|6.3.2| for a discussion of the systematic errors. The large 
symbols and navy blue plusses/crosses show the positions 
and metallicitie s for open clusters and Cepheids from the 
literature (see {7.2). 

In the low \Z\ slices {\Z\ < 0.5 kpc), we obtain values 
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TABLE 4 

Open Clusters with High-Resolution Observations. 



Cluster 


I 


b 


d (kpc) 


R (kpc) 


Z (kpc) 


[Fo/H] 


Age (Gyr) 


Bol7 


176.0 


-3.6 


2.7 


10.7 


-0.17 


-0.2 


10.1 


Be20 


203.5 


-17.3 


8.6 


15.9 


-2.55 


-0.5 


4.1 


Be20 


203.5 


-17.4 


8.3 


15.6 


-2.49 


-0.3 


6.0 


Be22 


199.9 


-8.1 


6.0 


13.7 


-0.84 


-0.3 


3.3 


Be25 


226.6 


-9.7 


11.4 


17.7 


-1.92 


-0.2 


5.0 


Bo29 


197.9 


8.0 


14.8 


22.4 


4-2.05 


-0.5 


4.3 


Be29 


198.0 


8.0 


13.2 


20.8 


4-1.83 


-0.3 


4.0 


Bc29 


198.0 


8.1 


13.2 


20.8 


-1-1.85 


-0.4 


4.5 


Be31 


206.2 


5.1 


5.3 


12.9 


-1-0.47 


-0.5 


5.3 


Be32 


208.0 


4.4 


3.4 


11.1 


-1-0.26 


-0.3 


5.9 


Be32 


208.0 


4.4 


3.1 


10.9 


-1-0.24 


-0.3 


5.5 


Be39 


223.5 


10.1 


4.3 


11.4 


4-0.75 


-0.2 


7.0 


Be73 


215.3 


-9.4 


9.8 


16.8 


-1.60 


-0.2 


1.5 


Be75 


234.3 
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th at are consistent wi th the —0.06 dex kpc ^ determined 
by Friel et al. ( 2002 1 for open clusters and iLuck et al. 



(2006) tor cepheids. Our major result is that in the high 
\Z\ slices {\Z\ > 0.5 kpc), the slope is flat for the entire 
radial range 6 < i? < 16 kpc. The constant median 
[Fe/H] in the highest slice {\Z\ > 1.0 kpc ) is consistent 



with the metallicities reported by Yong et al.| ( 2005 ) and 
others at large R. 

Our results are summarized in Figure |8j which shows 
the measured slopes as a function of vertical height \Z\. 
The horizontal dotted lines show the values in the lit- 



erature published by Friel et al. (2002) and Luck et al 



( 2006 ) for open clusters and cepheids, respectively, which 



are consistent with the values we obtained for the low \Z\ 
slices. Again, the red squares show the values obtained 
after we apply our weights, and the blue triangles show 
the values obtained using the unweighted data without 
any corrections for the selection function. The observed 
trend — flattening slope with increasing \Z\ — is seen in 
both the unweighted and weighted cases. 

To test the robustness of this result, we re- measured 
the gradient for various subsamples of our data. The 
purple symbols in Figure |8] show the results if we exclude 



likely halo stars (see ^6.2). The gray symbols show the 
values obtained for restricted ranges in R (7.3 < R < 
12.6 kpc and 7.8 < R < 10.5 kpc). These measurements 
were done to test that the flattening of the radial gradient 
is truly a trend in \Z\, and not just a result of the wider 
range in R that is observed in the high \Z\ bins. The 
two restricted ranges in R correspond to the extent of 
observations in the two lower \Z\ bins. Figure |8] shows 
that gradient becomes flat with or without weighting and 
regardless of the range in R observed. 

6. ERRORS 

6.1. Metallicity Bias 

As described in §2.1[ our sample of MSTO stars is 
chosen by making a cut in (g — r)sFD, which likely bi- 
ases it against redder, more metal-rich stars. Theoretical 
isochrones give us an idea of the significance of this bias. 
Figure^shows the temperature and (5 — ?')sfd distribu- 
tions oimetal-rich stars in our sample ([Fe/H] > -1- 0.2, 
gray circles ) against the theoretical isochrones of An 
et al.| ( |2009[ color). We estimate that the maximum ob 
servable ages of stars at [Fe/H] = -fO.4 and -1-0.2 are ~ i 
and 11 Gyr, respectively. At each metallicity, this corre 
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Fig. 7. — Metallicity [Fe/H] vs. Galactocentric radius R in four \Z\ slices. Light blue points indicate the SEGUE data. The weighted 
median metallicity and the derived linear fit are shown as red squares, with the numerical values in the bottom left of each panel. The 
blue triangles and values in parentheses show the results we would have obtained if no corrections for known selection effects had been 
applied. The spacing of the symbo ls ind icates the radial distribution of the targets. Open symbols and plusses/crosses are open clusters and 
Cepheids from the literature (see i ]7.2| for details). The sizes of the open cluster symbols indicate their ages (smaller symbols for younger 
clusters). At low \Z\ (< 0.5 kpc) our derived gradi ent is consist ent with published values. At high \Z\ (> 0.5 kpc) the constant [Fe/H] is 
consistent with the cluster metallicities reported by |Yong et al. |(2005) in the outer disk. 



spends to the turnoff age of the coolest /reddest stars that 
pass the selection criteria for our sample. The TAMT 
and maximum observable age are shown as solid lines 



Based on the measured ages o f nearby stars ( Bensby et al. 
2005 Casagrande et al.|2011 ), very few stars with metal 



licities [I'e/HJ > +U.2 are older than 8 Gyr. We therefore 
expect that metal-rich stars are well represented in our 
sample. For stars with metallicities [Fe/H] < 0.0, there 
is no significant bias against stars of any age. 

The arguments presented above suggest that the 
metallicity bias resulting from our color selection does 
not eliminate a significant fraction of the disk population 
from our sample. However, the color selection also causes 
metal-rich stars to be under-represented in our data. We 
use mock c atalogs from the models of [Schonrich fc Bin-| 
ney ( 2009a ) to quantify how well we correct for this effect 
and recover the true r netalli city distribution in our sam- 
ple. We show below {{ 6.3.2 ) that we are able to measure 
the true gradients within the errors and reproduce the 
trends as a function J_Z| when we apply the weighting 
scheme described in ^Mland the Appendix. 

6.2. Halo Contamination 



Since we are interested in the metallicity distribution 
of the Galactic disk, we must assess whether our sam- 
ple may be contaminated by halo objects. Our sample 
does not reach R < 5 kpc, where the bulge is expected 
to be a significant population. One way to quantify the 
amount of contamination is by examining different multi- 
component models for t he Galaxy, in particula r, those of 
Juri6 et al.| (|2008[) and |de Jong et "ah] ( |2010[ ). For the 



lines of sight in our sample, these two models predict 
total halo contaminations of ~ 2% and ~ 0.8%, respec- 
tively. In the most distant \Z\ bin, the predicted con- 
taminations increase to ~ 11% and ~ 5.6%, respectively. 
The difference in the predictions can be almost entirely 
attributed to the discrepant local den sities for the halo 
th at were derived (0.51 % by |Juric et al.,2008, and 0.17% 
b yjde Jong et al.||2010[ ). 

Both models, however, predict that the anti-center 
lines of sight, aimed toward the outer disk, will be more 
con taminated by halo stars. This is especially true in 
the Juric et al. (20081 model, which predicts that the 
two lines ot sight that reach the largest values of R will 
have 17 — 18% contamination (compared to 9 — 13% in 
the other lines of sight). If there is indeed more halo 
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Fig. 8. — Measured gradient A[Fc/H]/Ai? vs. vertical height 
\Z\. The horizontal black dotted lines indicate a flat gradient and 
the gradient measurement published by Friel ct al. (2002) using 
open clusters. The horizontal bl ue dotted Im es indicate the gradi- 
ent measurements published by Luck et al.i(2006) using Cephcids 
for the ranges 6.6 < i? < 10.6 kpc (zone ll) and 4.0 < i? < 14.6 
kpc (their total sample). The red squares and blue triangles show 
the measured gradients (weighted and unweighted, respectively) 
as a function of height above the Galactic plane presented in this 
work. The slope becomes flat at high \Z\ using both the weighted 
and unweighted results. The number of objects in each \Z\ slice 
is indicated above the symbols. The gray symbols show the re- 
sults obtained for restricted ranges in R. The purple symbols show 
the r esult s obtained when probable halo contaminants are removed 
(see i ]6.2| . In each case, the trend of the gradient becoming more 
shallow with distance from the plane is still observed, which indi- 
cates that our result does not arise due to the increase in the radial 
range or contamination from the halo. 

contamination in the outer regions, this should push our 
measured gradient to be steeper, since there will be more 
metal-poor halo stars at large R compared to the inner 
disk. However, we see the opposite effect, that the gra- 
dient is flatter in the highest \Z\ bins. 

Looking at the data, we do not see evidence for signif- 
icant halo contamination. Halo stars, which have large 
velocities and are more metal poor than disk stars, can 
generally be identified using their kinematics or chem- 
istry. To assess the effect that such stars may have on our 
results, we recalculate the gradient after applying two dif- 
ferent cuts to remove potential halo stars from the data. 
(1) A metallicity cut that removes all stars with [Fe/H] 
< —0.7, and (2) a kinematic cut that removes all stars 
with Vcai < 100 km s^^ to remove stars with the largest 
velocity offset relative to the projection of the local stan- 
dard of rest, where Vcai = Vu + 220 • cos6 • sinZ and Vr 
is the line-of-sight velocity measured from the SEGUE 
spectra. We only remove stars with Vcai < 100 km s^^ 
in lines of sight with 50 < Z < 130°. We do not include 
the lines of sight directed toward the Galactic anti-center 
because the local standard of rest is tangent to those di- 
rections and the projection does not give a meaningful 
velocity. The two cuts remove ~ 700 and ~ 100 stars, re- 
spectively. The gradients we measure using each of these 
cuts are not significantly different from our main results, 
as shown by the purple symbols in Figure [8] 

6.3. Mock Catalog Analysis 

To quantify the errors in our analysis, we utilize mock 
observa. t ions ge nerated from the models of |Schonrich &:| 
Binney (2009a I. By using a model where we know the 
true stellar parameters, distances, and metallicity dis- 
tributions of the targets, we can test whether we are 



able to reliably reproduce the ground truth after apply- 
ing our observational selection to the mock catalogs and 
performing the same analysis procedures. As the pur- 
pose is merely to assess how accurately we can measure 
quantities such as distances and metallicity gradients, 
this way of testing our methods is not dependent on 
having a correct model for the Galaxy. We do, how- 
ever, need a model that can reproduce som e basic prop- 
erties o f the ob served stellar populations. Schonrich fc 
Binney (2009a I have shown that their model provides a 
good niatch to the properties of stars in the solar neigh- 
borhood as observed by the Geneva- Copenhagen Survey 
( [Nordstrom et aL||2004[ ). 

We use a mock catalog, provided by R. Schonrich, of 
6,701,170 stars with ugriz colors, stellar parameters, and 
distances along ten lines of si ght. We repeat the same 
sample selection described in §2.1| on the mock catalog 
to replicate the effects of using MSTO stars as tracers. 
Figure [TO] is analogous to Figure [T] and shows the results 
of the random subsampling for two lines of sight. Though 
the model was not tuned to look like the SEGUE data, 
the model CMDs are good matches to the observations. 

In our analysis, we use all 111,640 objects that fulfill 
t he MS TO selection to estimate errors on the distance 
(^6.3.1). We draw many random subsamples of 6500 
MS'l'O stars to estimate the systema tic an d random er- 
rors on our gradient measurement ([6.3.21, to simulate 
the effect of having a limited number of spectroscopic 
targets in each line of sight. We repeat the analysis us- 
ing different random subsamples to assess how much the 
results change based on which particular targets are cho- 
sen for spectroscopy. 

6.3.1. Errors in Distance Estimates 

Systematic Errors: The main source of systematic 
error in the distance arises from assumptions we must 
make when choosing isochrones to estimate the lumi- 
nosities of o ur stars. One s ource of error is the a- 
enhancement; An et al. (2009j) adopt an a-enhancement 
scheme where each value ol [i'e/H] has an assumed value 
of [a/Fe] (see Table 2). We do not expect the discrep- 
ancy between the target and isochrone [a/Fe] to have a 
large effect on the distance estimate. The [a/Fe] of stars 
in our sample ( jLee et al.| |2011b[ ) are generally within 
0.2 dex of the values assigned to the An et all (2009[) 
isochr ones. Using the [Fe/H]-[a/Fe] grid of Uotter et al. 
( |2008f we estimate that a 0.2 dex change in [a/Fe], at 
worstjTeads to a ~ 10% change in the distance for a star 
on the zero-age main sequence. 

For our sample, a larger source of uncertainty is the 
ages we assign to our targets; the slope of the main se- 
quence becomes very steep at the turnoff, making the 
predicted distance a strong function of the chosen age. 
We compare the distance estimates obtained using three 
different age assumptions, which are shown schematically 
for stars at solar metallicity (—0.05 < [Fe/H] < 0.05) in 
Figure [3j 

1. Turnoff Age of the Mean Temperature (TAMT): 
We use one isochrone for most of the stars in the 
metallicity bin — this is determined by finding the 

^ http:/ /stellar. dartmouth.edu/~models/grid. html 
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Fig. 9. — Maximum observable ages for metal-rich stars. Wo compare the SEGUE data (gray circles) with theoretical isochrones (color, 

|An et al. 2009 ) to estimate the oldest main sequence turnofT stars that would fall in our sample. The TAMT (see Table Hi and maximum 
oBservable age are shown as solid lines. For [Fe/H]= +0.4 (left panel), the coolest /reddest stars may be as old as ~ 8 Gyr; lor [Fe/H]= +0.2 

i right panel), the coolest /reddest stars that fall in our sample may be as old as ~ 11 Gyr. Based on the measured ages of nearby stars 
Bensby et al.|2005[|Casagrande et al.|2011| , we do not expect there to be a significant population of old, metal-rich stars which would be 
excluded trom our sample. 
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Fig. 1 0. — Mock catalog MSTO selection. This figure is analogous to Figure[T] but using simulated observations from |Schonrich &: Binney] 
(2009aJ. The full sample of stars (grayscale and contours) is put through the same sample selection as the real data to yield an analogous 
Mb'l U sample (blue circles). The similarity between the CMDs of the real and mock data shows that we are using a reasonable model of 
the Galaxy to test our errors. 



mean of the temperature distribution. Targets hot- 
ter than the main sequence turnoff are assigned to 
the oldest possible isochrone. This was the scheme 
used throughout our analysis (described in fj3| and 
illustrated in the second panel of Figure |3] 

Turnoff (TO): We use the oldest possible isochrone 
consistent with their measured temperature for all 
targets — this assumes that every star is located at 
the turnoff. In comparison to the TAMT assump- 
tion, this approach changes the distances for stars 
cooler than the mean temperature of the sample, 
as shown in the third panel of Figure |3j 

Zero- Age Main Sequence (ZAMS): We use the 
youngest possible isochrone for all targets — this as- 
sumes that most targets are located on the zero-age 
main sequence, as shown in the fourth panel of Fig- 
ure |3l 



For this analysis, we compare our calculated distance 
d to the true distance from the model dmodei for 111,640 
mock targets at solar metallicity (—0.05 < [Fe/H] < 
0.05), where theoretical isochrones from different groups 
show the best agreement. Any disagreement between the 
isochrones used to generate the mock catalog and the 
isochrones used to calculate the distances will introduce 
an additional systematic error. The distributions of the 
systematic fractional error {Ad/dhys for all three age as- 
sumptions are shown in Figure |11[ where the systematic 
fractional error is given by: 



Ad 
~d 



d-d 



model 



sys 



odd 



(1) 



Based on the offsets of the distributions, the TO (red 
dashed line) and ZAMS (blue dash-dotted line) assump- 
tions overestimate and underestimate the distances, re- 
spectively, though the ZAMS assumption does much 
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Fig. 11. — The fractional error distributions of solar metallic- 
ity stars (—0.05 < [Fe/H] < 0.05) for three different age assump- 
tions, as described in the text (see also Figure [sjl. The TAMT 
age assumption gives the best agreement between the calculated 
and true model distances. Based on the width of the distribution, 
we estimate the systematic error on our distance estimates to be 
~ 10%. The vertical line indicates where the fractional error is 
zero. The gray line is a Gaussian fit to the TAMT distribution 
where —0.15 < (A(i/(i)sys < 0.0, which we use to estimate the 
fraction of subgiants in our sample (~ 15%). 



worse. The TAMT (black solid line) assumption gives 
the best agreement between the calculated and model 
distances, with an error of ~ 10%, based on the width of 
the (A(i/(i)sys distribution. Mock targets at other metal- 
licities exhibit the same behavior, with the TAMT as- 
sumption giving the best distances with a mean error 
of ^ 10%. This error includes the discrepancy between 
the isochrone [a/Fe] and that of individual stars, as the 
model stars have a range of [a/Fe] that d o not exactly 
match the sequence of the An et al. (2009) isochrones. 

Using the (Ad/d)sys distribution, we can also estimate 
the number of subgiants in our sample by investigating 
the lon g ta il of negative fractional errors. The gray line in 
Figure 11 shows a Gaussian fit to the TAMT histogram 



where —0.15 < (Ac?/(i)sys < 0.0 — this range was chosen 
by eye to reflect the range where subgiants were not con- 
tributing to the distribution. By examining the discrep- 
ancy between the Gaussian fit and the long tail of nega- 
tive (Aci/(i)sys values, we estimate the contamination by 
subgiants to be ^ 15%; we obtain the same result if we 



sys 



< 



-0.15. 



take subgiants to be all stars with (Ad/c?) 
We derive similar values for other metallicities as well. 

We can also assess how an error in the value of the 
TAMT affects the distances that we obtain in our cal- 
culation. If we view the TO assumption as the TAMT 
assumption with an incorrect TAMT — one that is too 
old by 0.2 dex — we can estimate how much the distances 
change as a result. Comparing the TAMT and TO re- 
sults in Figure 11 the peak (Ad/c?)sys values of the two 



assumptions are within 5%, which suggests that any er- 
ror in the TAMT is smaller than the error from the age 
assumption scheme used. 

Random Errors: The random error in the distances 
is dominated by uncertainties in the stellar parameters 
derived by the SSPP; errors in the magnitudes are triv- 
ial in comparison. We examine how the parameter errors 
propagate through our analysis by generating 500 Monte 
Carlo realizations of the mock catalog with slightly per- 



turbed values of the stellar parameters and repeating the 
same distance measurements each time. For each realiza- 
tion we assign a Gaussian distribution of perturbations 
with widths of 200 K and 0.3 dex for Tog and [Fe/H], 
respectively. These values are motivated by the errors 
estimated by Lee et al. (2008a). The resulting random 
fractional error (Ad/d)rand, given by: 



Ao?\ _ dpcrtuibcd 

^ Aand ^ 



(2) 



where dporturbcd is the distance calculated with the new 
perturbed values of the stellar parameters. From the dis- 
tribution of (Ad/(i)i.and, we estimate that the uncertainty 
in the SSPP parameters translates to a ~ 15 — 20% error 
in the distance. These values correspond to the 68% con- 
fidence levels derived from the width of the distribution. 

We can now calculate the total error in the distance, 
which combines the systematic and random errors from 
the age assumption and the stellar parameters, respec- 
tively. The result is the total fractional error (Ad/d)totai, 
given by: 



Ad 



d, 



perturbed 



^model 



total 



^modcl 



(3) 



We estimate that including both the systematic and ran- 
dom errors, we have a total distance error of ~ 20—25%, 
with a larger contribution from the uncertainty in the 
stellar parameters determined by the SSPP. This total 
error is comparabl e to the size of the e rrors in the dis- 
tances estimated by Gilmore et al. ( 1995 ) for their sample 
of MSTO stars. 

6.3.2. Errors m the Metallicity Gradient Measurement 

Systematic Errors: Using the same mock catalog, 
we also assess how well we measure the radial metal- 
licity gradients in the disk. Figure 12 is the same as 
Figure [8] but shows the gradient resultslor the mock cat- 
alogs in four |Z| slices. The true gradients, measured 
using the entire mock catalog, are shown as black circles. 
We draw 100 random samples (6500 MSTO stars each) 
from the mock ca talo g using the same selection criteria 
as the real data {\ 2.1 ) and follow the same analysis pro- 
cedures, including accounting for the weights (Mland the 
Appendix). 

The filled symbols in Figure [12] show the mean gradi- 
ent of the 100 realizations; the error bars indicate the 
standard deviation of the distribution. The smaller open 
symbols show the gradient measurements for one partic- 
ular realization, with the error bar indicating the random 
error due to the uncertainty in the stellar parameters (see 
below) . 

The true gradients (open black circles) are generally 
within the errors of the weighted measurements (filled 
red squares), indicating that our measurements are re- 
liable. More importantly, the weighted measurements 
show the same flattening trend with |Z| as the true gra- 
dient. This is not seen in the unweighted gradient mea- 
surements (filled blue triangles). While the metallicity 
gradient in t he model is steep compa red to recent obser- 



vations (e.g.. Luck & Lambert 2011), our mock catalog 



analysis demonstrates that the weighting procedure is 
necessary to reproduce the attening trend. 
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Fig. 12. — Measured gradient A[Fe/H]/A_R vs. vertical height 
\Z\ using mock samples. This figure is analogous to Figure [s] ex- 
cept the samples used are dra wn from mock catalogs of the galaxy 
model of [Schonrich Binney| ( j2009a^). The mean and standard 
deviation tor iUU sucii samples ot Mb t O stars are shown as filled 
symbols. The open symbols show the measured gradients for one 
particular random sample, with the error bars showing the random 
errors derived from Monte Carlo realizations using perturbed stel- 
lar parameters. The black circles show the true gradient, which is 
measured using all of the targets in the catalog. The true gradient 
generally falls within the errors presented for the weighted gradient 
(red squares) measured using the MSTO samples, which indicates 
that our method of measuring the gradient gives a reliable result. 
More importantly, we are able to reproduce the flattening trend 
seen in the true gradient. This is not true for the unweighted gra- 
dient measurements (blue triaiigles). Note that the vertical scale 
is different from that of Figure [8| 

Random Errors: The error bars on the open symbols 
represent the random errors, which are derived by run- 
ning our full analysis on one sample of 6500 MSTO stars 
from the mock catalog, where we generate 500 Monte 
Carlo realizations of the mock data with perturbed stel- 
lar parameters. The error bars indicate the 68% confi- 
dence levels derived from the width of the distribution of 
slopes we obtain in the 500 realizations. This is the same 
procedu re foll owed to determine the random distance er- 
rors in §6.3. 1| The random errors are comparable to the 
systematic errors on the gradient, although in some cases 
they are larger. The large error bars for the lowest \Z\ 
bin is due to a combination of the smaller number of stars 
and the smaller range in R. 

7. DISCUSSION 

7.1. Comparison with Previous Studies 

The median metallicity we find at |Z| > 1 kpc 
([F e/H] ^ —0.5) is cons istent with the value published 



by Gilmore et al. (1995) for their stars located 1.5 kpc 



above the plane; the median \Z\ for our sample above 1 
kpc is 1.24 kpc. Their sample also consisted of an in-situ 
sample of F/G dwarfs, a llowing for a direct comparison. 

also drawn 



J 



2011a 



The G dwarf sample of Lee et al 
from SEGUE, has median values 
dex of that of our sample as a function of 

Our results are also consistent with the lack of a ra- 
dial m etallicity gradient found by both 



that we observe lower Galactic latitudes and we are able 
to make a direct comparison of the thin and thick disks 
using the same homogeneous sample. Our results are 
based on a sample that has been carefully corrected for 
selection effects and we use an improved version of the 
SSPP that is more accurate for metal-rich stars, which 
was not available for these previous analyses. 
In their analysis of thick disk stars i n the solar neigh- 



2005[ ) find that abun- 
and b'e do not vary as a 



borhood, Bensby et al. (2003 
dance trends using U, Mg 

function of R or l^maxl in the region 5 < R < 7 kpc, 
< l^maxl < 1-1 kpc, where |^max| is the maximum 
vertical height reached by the calculated stellar orbit. 
These detailed abundances, together with our finding of 
a flat metallicity gradient at large \Z\, are suggestive of 
a chemically homogeneous thick disk, although more ob- 
servational data are needed to confirm this. 

7.2. Comparison with Open Clusters and Cepheids 

In Figures [5] and [7] we show open cluster and Gepheid 
data from the literature. These two classes of objects 
have been studied extensively with high-resolution spec- 
troscopy and also span a large range in age, giving them 
the power to probe the chemistry of the interstellar 
medium at different times. Cepheids typically have life- 
times on the order of ~ 100 Myr and trace the present 
day metallicity distribution. Open clusters can have ages 
anywhere from less than a Gyr to greater than 10 Gyr 
and have been used to examine the time evolution of the 
disk's metallicity gradient. 

We compare our old disk stars to 33 open clusters with 
abundance determinations from high-resolution spectra. 
These data are obtained from the work of four groups 
that h ave studied clusters at i? > 10 kpc : 5 from|Yong 



et al. (12005, Y05), 9 from Carraro et al. (Carraro et al. 



^nn^lMFTriVillanova et al. 2005; Frincha boyet al.|i0081 



C07), 1 2 from Bragagli a, Carretta, Sestito7et al 



gagha e t al. 2001; Carre tta et al.|2004[ [20051 [2007| |S]estito 

13 from Friei, Jacob 



et al.||2006| 



son, and 



6i |2007n2008| BGS), and 
Piiacho wski | Fn el et al. 



2005 2010 Jacobson 



et al.| |2008; '2009V2011"I''JP). Seven clusters have mea- 
sured abundances Irom more than one study, which pro- 
vides an indication of the size of the uncertainties and 
systematics between groups. Properties of the clusters 
are listed in Table ^ The open clusters have ages rang- 
ing from 0.6 to 10.1 Gyr; the median age is 4.0 Gyr, 
with only two as old as 7 Gyr, which makes this clus- 
ter sample younger than the typical age of the stars in 
our sample (see Table [s]). The ages of the clusters are 
indicated by the size oithe symbols in Figure [5] and [Tj 
with older clusters having larger symbols. For clusters 
with more than one measurement, a line connects the 

symbols representing the different det erminations . 

The Cepheid d ata are taken from [Andrievsky et al 



JFejk] within 0.1 (2002c a b, 2004] ; [Hlcketal] ([2003 



|et al.| ( j2006t Y05). Most of the 



2006 



et al. 



Allende Prieto 



( |2006J , using spe c tra of 12,483 F/G stars from the 



SPSS, and Juric et al.[ ( [2008[ ), using photometric metal- 
licities for millions ot stars m the SDSS, but both stud- 
ies examine a much higher \Z\ sample with a limited R 
range at the vertical heights that our sample covers. The 
present study is complementary to these two studies in 



AL), and Tong 

lieids In these samples 
are too close to the midplane {\Z\ < 0.15 kpc) to be 
directly compared to our sample. Those in the higher 
bins (0.5 < \Z\ < 1.0 kpc), which are mostly from the 



Yong et al. (2006) sample, tend to be at larger radii than 
our MS'l'U stars. For consistency we have re-calculated 
R, \Z\ values using published I, b, and distances with 
Rgc,q = 8.0 kpc where necessary. 
Both Cepheids and open clusters appear to be slightly 
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more metal-rich than the median metaUicities of the old 
disk stars. But systematic differences in [Fe/H] between 
different groups can be up to 0.2 dex, as shown by the 
clusters with multiple measurements ( see Table p] as wel l 
as the discussion in the Appendix of Friel et &l\ 2010). 
We do not see any obvious differences in metaUicities be- 
tween old and young clusters (large versus small open 
symbols). At low \Z\, the slopes of the open clusters and 
field stars do not appear to be dramatically different, as 
expected given the good agreement between our measure- 
ment of the radial gradien t in the low \Z\ bins and the 
values o flFriel et al.| ( |2002[ ) and |Luck et al.[ ( |2006| dotted 
lines in Figure [8]). However, at high \Z\, there are three 
clusters with metaUicities ~ 0.5 dex higher than the me- 
dian [Fe/H] for field stars at small R (< 10 kpc). While a 
steeper slope in the inner regions of the disk (—0.13 dex 
kpc~^ at i? < 8 kpc) has been reported for both open 



to large R (> 15kpc) for all \Z\, then the discontinu- 
ity in the radial metallicity gradient of the d isk actually 



reflect s the change in slope at different \Z\. Friel et al. 
(20101 have stressed the need for more high-resolution 
data of open clusters to build up a homogeneous, statis- 
tical sample; our work shows that it is also necessary to 
fully sample a range of both R and \Z\ to understand 
the metallicity trends in the disk. The large R, low \Z\ 
region of the Gala xy will be probed by future surveys 
such as APOGEE ( [Eisenstein et al.||2011[ ). 



clusters and Cepheids (e.g., Pedicelh et al.j200l) Magrini 
eraI1[20T0| , we do not see a steeper metallicity gradient 



7.3. Implications for Thick Disk Formation 

The lack of a radial metallicity gradient far from the 
Galactic plane provides an observational constraint that 
must be matched by any viable scenario for thick disk 
formation, such as the four described in Sjl] A flat gra- 
dient, as we measure, may come most naturally out of a 
turbulent disk at high redshift (scenar io 3). If the thick 



at R < H kpc in our sample of old disk stars 

The question of how the radial metallicity gradient may 
change over time is still debated in t he literature , with 
conflic ting reports of fl attening (e.g. 
Chen et al. 2003 



disk f ormed rapidly at early times (e.g., Brook et al. 2004 
2005 1 , then the thick disk would be chemically homoge 



Friel et al.||2002 

, Mac iel et al. 2005), steepenmg (e.g., 

Stanghellini fc Ha y wood | 2010 ), and constant slopes (e.g., 
Magrini et al. 2009| as a tunction of time. If the gradi- 
ent IS steeper m the inner disk, as has been reported, the 
flatness of the gradient for old disk stars, at face value, 
implies that the radial metallicity gradient in the inner 
disk has grown steeper with time. Radial migration com- 
plicates this issue greatly, however, as the movement of 
stars from their initial orbits may wash out a previously 



neous and the radial metallicity gradient far from the 
plane of the disk would be flat. 
The observations can only be explained by current 



models of radial migration in isolated di sks (e.g., Roskar 
et al.|[20b8a Schonrich fc Binney|2009a I if radial mixing 
is strong (scenario 4). 'I'hese models show that dynam- 
ical interactions with spiral arms can move stars from 
their initial orbits and make the gradient shallow, but the 
mechanism must be efflcient if it is the dominant mech- 
anism for forming the thick disk. The degree of radial 
mixing may be increas ed through dynamical interactions 
with the Galactic bar ([Friedh e t a L]|1994[ [Martin fc Roy 



the observations presented here, we are not able to draw 
any clear conclusions about the time evolution of the 
metallicity gradient in the inner disk. 

The disagreement in abundances between young and 
old tracers is not present in the outer disk, where the 
median metallicity of the old disk stars ([Fe/ H] ~ —0.5) 
is consistent with the metaUicities reported by |Yong et al.| 
([2005) for outer disk open clusters. Furthermore, the me- 
dian metallicity of our sample at \Z\ > 1.0 kpc is constant 
at all values of R, which suggests that the flattening of 
the gradient is due to a trend in \Z\. Given the change in 
the radial range spanned in each \Z\ bin, a superposition 
of the negative gradient at small \Z\ with a flat gradient 
at large \Z\ could result in an apparent discontinuity in 
the radial gradient. 

This result suggests that it is important to distinguish 
between trends in the radial and vertical directions. As 
shown in Figure [5] all of the outer disk clusters from 
the literature (i? > 15 kpc) are located far from the 
Galactic plane, and their mean metallicity is consistent 
with our field stars located at similar vertical heights 
{\Z\ > 1.0 kpc), although our field star sample extends 
only to i? ~ 15 kpc, whereas the clusters reach radial 
distances of 20 kpc or more. Based on our observations of 
old disk stars, whether the reported discontinuity is truly 
a feature of the radial metallicity gradient is unclear. 

Our sample does not reach to sufficiently large radii 
at low latitudes; with the limited range in R at low \Z\, 
we cannot confirm that the trend is purely vertical. If 
the radial gradients we measure at \Z\ < 1.0 kpc extend 



existing gradient (e.g., |Roskar et al.||2008a[). Based on |1994| [Minchev fc Famaey||2010p . 



Another possible way to increase the amount of radial 
mixing in a disk is to place the disk within a cosmologi- 
cally motiv ated accreti on history (scenario 1). The sim- 
ulations of Bird et al. ( 2011) show that stars can move 
far from their initial orbits in a disk that is bombarde d 
by multiple minor mergers (e.g., Kazantzidis et al.|2008 ). 
Furthermore, the bombarded disk experiences more ra- 
dial mixing at high |Z| than an identical disk in isolation. 
At the midplane, however, both scenarios show about the 
same amount of mixing. The increased radial mixing at 
high |Z| could be responsible for the lack of a radial gra- 
dient. 

With our data alone, we cannot rule out the direct ac- 
cretion of stars that originally formed in satellites that 
have been disrupted (scenario 2). To assess whether the 
thick disk exhibits the predicted dumpiness, we would 
need to investigate the azimuthal variation within our 
sample or at more detailed abundances with foUowup 
observations. But observati ons of the orbit al properties 



of stars in S PSS and RAVE ( jPierickx et al.l2010[|Wils(m| 
et al.|[20TT ) show that there are not enough higli eccen- 
tri city stars to match what is seen in the simulations 
oflAbadi et al.'('2003 ). Furtherm ore, recent simulations 
(Read et al. 2008: V ilfalobos fcHe lmi 2008) predict lower 
numbers of stars being directly accreted during minor 

mergers than the 2003 simulations. 

Our c onclusions are consistent with those of ILee et al.l 
(2011a), who favor a cosmological formation scenario for 
the thick disk; they base this interpretation on correla- 
tions between rotational velocity and [Fe/H] and R. 
Only recently have simulations begun to include prescrip- 
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tions for the formation of stars and their impact on the 
expected chemistry of the disk populations (e.g., Brook 

— ■ — ■ ■ — mi 



et al.|[2005l [Stinson et al. 



1 [20T0l 

he mec 



Loebman et al 



and references therein). The measurement of the radia' 
metalhcity gradient at different heights above the plane 
is an additional observational constraint that can be used 
to test new, improved models for disk evolution and the 
next generation of simulations. 



8. SUMMARY 

Using a sample of 7010 MSTO stars from the SEGUE 
survey, we measure the radial metalhcity gradient 
A[Fe/H]/Ai?, as a function of height above the plane \Z\, 
in the Milky Way disk. Near the midplane, where our 
sample is dominated by the thin disk, we see a negative 
radial metalhcity gradi ent that is consistent with previ 



radial metaliicity gradi ent tnat is consistent witn previ- 
ously published values ( Friel et al.|2002 Luck et al.|2006 



Luck & Lambert 



2011 1. At large vertical heights, where 



we are dominated by the thick disk, the radial metalii city 
gradient becomes flat, c onsistent with previous work ( Al 



lende Prieto et al. 12006) using a sample located at larger 



\Z\. Our sample, located at low Galactic latitude, covers 
a larger range in R at small \Z\ and allows us to make a 
direct comparison between the thin and thick disks using 
the same sample. 

At \Z\ > 1.0 kpc, the median metaliicity of old disk 
stars is co nsistent with the o pen cluster metallicities re- 
(2005[) and others at large R. In 
stars shows that the flat 



ported by Yong et al 
addition, our sample ot dis 
gradient at large vertical height \Z\ extends to small R. 
Because the outer disk clusters are all located at large 
\Z\, the reported discontinuity in the radial gradient is 
consistent with a transition found using tracers at small 
\Z\ to large \Z\. We stress that abundances need to be 
examined as a function of both R and \Z\ in order to 
truly understand the observed trends. 

In contrast to the outer disk, open clusters and 
Cepheids in the inner disk at high \Z\ have median metal- 
licities ~ 0.5 dex higher than old disk stars at the same 
R and \Z\; thus, far from the Galactic plane, the younger 
tracers do not exhibit the same flat metaliicity gradient 
that is seen in the old disk stars. Whether this is indica- 
tive of a metaliicity gradient that becomes steeper with 
time is unclear, as radial migration may play a role in 
erasing a pre-existing gradient in the old disk stars. 

A flat radial metaliicity gradient in the thick disk is 
consistent with the predictions of a gas-rich, turbulent 
disk at high redshift ( Brook et al.||2005 Bournaud et al. 
20091. It may al so be consistent with the scenarios of 
radial migration (IRoskar et al. "2008a Schonrich & Bin- 

tr 2009b; Minchcv & I'amacy 201( j| or minor mergers 
azantzidis et aTl^ODSj |Bird et al.||2t)ll[ ) , provided that 
mixing in the radial direction is strong. We also cannot 
exclude the direct accretion of stars from satellites in mi- 
nor mergers ( |Abadi et al.|2003| ). While we are not able to 
conclusively rule out any ot tnese scenarios, the change 
in the radial gradient as a function of height above the 
plane is an important observational constraint for future 



theoretical work. 
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APPENDIX 
APPENDIX: CALCULATING WEIGHTS 

As discussed in S|4j we assign weights to each of the MSTO stars to reconstruct the properties of the underlying 
parent population. There are three major ways in which the spectroscopic sample is different from the full population 
along each line of sight: (1) The photometric objects in regions with the highest extinction were not considered for 
spectroscopy. (2) Not all candidates for spectroscopy are observed. (3) We observe only MSTO stars using a color cut 
which is biased against redder metal-rich stars. 

Each star in our sample is given three weights which correspond with the three differences listed above: (1) the 
area weight Wa, which depends on the coverage on the plane of the sky in each line of sight; (2) the CMD weight 
WcMD, which depends on the target's location in the CMD; (3) the LF weight Wlf, which depends on the target's 
Toff, [Fe/H], and location in the CMD. The total weight W is the product of the three weights Wa, Wcmb, and Wlf- 

Area Weight 

The area weight Wa corrects for the area not covered by the spectroscopic survey because the 25% most extincted 
photometric objects were not considered for spectroscopy. These objects were removed using the cut described in steps 
2-3 of §2.1| The righthand panel of Figure |4] shows an example of how the missing area is distributed in the field. 
Since the missing area in each line of sight is slightly different, this weight is needed to ensure that every line of sight 
effectively probes the same volume of the Galaxy and has equal influence on the final gradient measurement. Since the 
dust is mainly in the foreground and our sample is primarily distant stars, we assume that the volume of the Galaxy 
behind the high extinction patc hes is the same as the rest of the volume probed by stars along the same line of sight. 

We use the extinction map of |SFD98| to calculate the area with extinction lower than the cut. Wa i s th e ratio of 
the total area (7 square degrees) to the low extinction area. We note that the angular resolution of the | SFD98, maps 
is 6.1 arcminutes, so that the most extincted regions are always in irregular, contiguous patches on the sky. 

Color- Magnitude Weight 

The CMD weight Wcmu normalizes between the different lines of sight; while each field has roughly the same number 
of spectra, the total number of photometric objects varies a great deal due to the structure of the Galaxy. This weight 
also accounts for any uneven sampling due to the stochastic nature of the random selection of spectroscopic targets, 
especially at the faint and red limits, where targets are less likely to end up in our samp le b ecause they have low 
qua lity spec tra. We divide the CMD into bins of gsFU and {g — r)sFD, as shown in Figure [l3| (black hues). We use 
the SFD98 colors and magnitudes because this is the CMD in which the g — r cut was applied, meaning that the 
ran domly se lected spectra are an unbiased sample of the underlying CMD; this procedure does not depend on whether 
the |SFD98| colors and magnitudes are correct. Since the sampling is a smooth function of color and magnitude and 
changes slowly, the bins are sufficiently small that we can assume the sampling is constant within each bin. 

The magnitud e bins ar e 0.5 magnitudes wide and span the entire range of the sample. Because the CMD bins use 
corrections from SFD98 which can be affected by reddening, the color bins are different for each line of sight. To 
determine the color bins, we calculate the mean and standard deviation of the (g — r)sFD colors of stars in each line 
of sight. Spectra with colors more than 2-a from the mean have VFcmd = 0; this removes 217 stars from the sample 
(open circles). Most of these are very blue objects (hotter stars that are likely not on the main sequence), although 
some red objects are removed as well. The remaining sample in the line of sight is divided into five equal-sized color 
bins. In each CMD bin, M^cmd is the ratio of the number of photometric objects (small gray dots) to the number of 
spectroscopic objects (blue circles) and is shown by the color coding in Figure |T3j The highest values of Wcmd are 
found in the reddest and faintest bins, while the variation in the middle of CMT) is relatively small. The difference 
between the smallest and largest values is less than a factor of ten. 

Luminosity Function Weight 

The LF weight Wlf allows us to use the MSTO sample as tracers for the total underlying population. The CMD 
bins used to calculate Wlf use the isochrone extinction corrected g^ and {g — r)o (see §3| and are different than 
those used to cal culate V Fcmd- For Wqmu, we had to account for the random sampling from the total photometric 
sample using the |SFD98| corrections, but for Wlf, we use the isochrone extinction corrections because they provide 
the best estimates of the absolute magnitude and g — r colors. The color bins are 0.1 magnitudes wide in the range 
0.05 < {g ~ f )o < 0.85, while the magnitude bins are 0.5 magnitudes wide in the range 12 < 50 < 20; these bins apply 
for all 11 lines of sight. 

In each CMD bin, we find the fraction of the luminosity function that is observed in th e given (q — r)n range. The 



weight is simply the reciprocal of the fraction. We use luminosity functio ns assuming a Chabrier (2001) lognormal 



initial mass function, generated by the Padova groupp] ( [Girardi et al. 2004), where we have modified the faint end of 
the theoreti cal luminosi t y fun ction s to more closely reflect the shapes ot the luminosity functions for disk field stars 
reported by 
the peak of 1 



Reid et al. (2002) and Bochanski et al.| (|2010 ). The modifications are made at magnitudes fainter than 
he luminosity function. 'I'he luminosity function used depends on the age and metallicity of the target. 
The derived values for Wlf are shown in Figure [14] for the entire sample; each panel shows a different range in 
metallicity. For a given {g — r)o color, Wlf is roughly constant as a function of magnitude. Within a given color bin. 



http:/ / stev.oapd.inaf.it / cmd| 
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Fig. 13. — Weighting for CMD sampling along one line of sight. For each bin (black lines) the CMD weight Wcmd is calculated by taking 
the ratio of the number of photometric (gray) to spectroscopic (blue) objects. Outliers bluer or redder than 2-cr from the mean value of 
(9 ~ '')sFD are given a weight of zero (open circles). The values of Wqmd vary by less than a factor of ten, with the highest values in the 
reddest and faintest bins. This correction accounts for uneven sampling of the CMD and provides the normalization for the variation in 
the number of objects between different lines of sight. 

the change in Wlf corresponds to different main sequence turnoffs for isochrones of different ages, which are assigned 
based on the temperatures of the targets as described in S|3] and Figure |3] In a given CMD bin, stars with older ages 
have larger LF weights because a smaller fraction of the luminosity function is observable within the bin. This accounts 
for the vertical striping pattern seen in the four panels. We assign Wlf = to the three spectra that fall outside the 
bounds of the CMD bins. 

Using the luminosity functions to correct for the fraction of unobserved stars represented by the stars in each CMD 
bin corrects for the bias against selecting metal-rich stars for the spectroscopic sample. The identification of MSTO 
stars relies on a cut in (g — r)gFD, which preferentially removes metal-rich stars that fall on isochrones with redder 
turnoff colors. Thus, for metal-rich stars, more of the population is removed by the color cut and a smaller fraction of 
the luminosity function is observed. Consequently nietal-rich targets have larger values of Wlf (lighter blue circles) 



at the same (17 — r)o in the four panels in Figure 14 In order for VFlf to make the proper corrections, we need to 
have a sufficient number of metal-rich stars in our sample to which we can apply the weights. Based on the arguments 
presented in |6.1[ our sample should satisfy this requirement. 



REFERENCES 



Abadi, M. C, Navarro, J. F., Steinmetz, M., & Eke, V. R. 2003, 

ApJ, 597, 21 
Abazajian, K. N., et al. 2009, ApJS, 182, 543 
Aihara, H., et al. 2011, ApJS, 193, 29 

AUende Prieto, C, Beers, T. C, Wilhelm, R., Newberg, H. J., 
Rockosi, C. M., Yanny, B., & Lee, Y. S. 2006, ApJ, 636, 804 
An, D., et al. 2009, ApJ, 700, 523 

Andrievsky, S. M., Bersier, D., Kovtyukh, V. V., Luck, R. E., 
Maciel, W. J., Lepine, J. R. D., & Beletsky, Y. V. 2002a, A&A, 
384, 140 



Andrievsky, S. M., Kovtyukh, V. V., Luck, R. E., Lepine, 
J. R. D., Maciel, W. J., & Beletsky, Y. V. 2002b, A&A, 392, 
491 

Andrievsky, S. M., Luck, R. E., Martin, P., & Lepine, J. R. D. 

2004, A&A, 413, 159 
Andrievsky, S. M., et al. 2002c, A&A, 381, 32 
Belokurov, V., et al. 2007, ApJ, 658, 337 

Bensby, T., Feltzing, S., & Lundstrom, I. 2003, A&A, 410, 527 
— . 2004, A&A, 421, 969 

Bensby, T., Feltzing, S., Lundstrom, I., & Ilyin, I. 2005, A&A, 
433, 185 



Metallicity Gradients Observed by SEGUE 



21 



12 
14 

1^ Ifi 

1& 

20 

12 
14 

^ ie 
ia 

20 









1 













































^ ■ 














** 




■Jm 

■ 


































































n-; 
























1 


L J J . -.1 


L * 




i , 



0,2 QA 0.6 0.B 




2.5 



2J ■ 



firi ■ .1 ■ ■ 

-'.5 1 

■HI.: '-'i - . 



* ■ I L 



12 
14 

A ie 

ia 

20 



0.2 0,4 0,fi 0.B 

































r 












■*p — 










- 


— 






t . 












1 


m 












■ 


.* r 




















■ 2 




* 
































































■1 
















¥— 






















■ 


1 










, 1 














0.4 




OA 




a5 



2.0 



1J 



1-S 



1.2 



T,Ol 



Fig. 14. — Weighting for LF coverage for tlie total sample. In each isochrone extinction-corrected CMD bin (black lines) the LF weight 
Wlf corrects for the fraction of the luminosity function that is observed. The luminosity function used varies with both metallicity (shown 
in four separate panels) and age, which depends on the temperature. This correction accounts for the bias against redder, metal-rich stars 
that results from the (g — »')sfd color selection. 



Bird, J. C, Kazantzidis, S., & Weinberg, D. H. 2011, ArXiv 
e-prints 

Bochanski, J. J., Hawley, S. L., Covey, K. R., West, A. A., Reid, 

I. N., GoUmowski, D. A., & Ivezic, Z. 2010, AJ, 139, 2679 
Bournaud, F., Elmegreen, B. G., & Martig, M. 2009, ApJ, 707, LI 
Bragaglia, A., Sestito, P., Villanova, S., Carretta, E., Randich, S., 

& Tosi, M. 2008, A&A, 480, 79 
Bragaglia, A., et al. 2001, AJ, 121, 327 
Brook, C. B., Gibson, B. K., Martel, H., 

ApJ, 630, 298 
Brook, C. B., Kawata, D., Gibson, B. K, 

ApJ, 612, 894 
Brooks, A. M., Governato, F., Quinn, T 

Wadsley, J. 2009, ApJ, 694, 396 
Caputo, F., Marconi, M., Musella, I., & 

544 

Carraro, G., Bresolin, F., Villanova, S., Matteucci, F., Patat, F., 

& Romaniello, M. 2004, AJ, 128, 1676 
Carraro, G., Geisler, D., Villanova, S., Frinchaboy, P. M., & 

Majewski, S. R. 2007, A&A, 476, 217 
Carretta, E., Bragaglia, A., & Gratton, R. G. 2007, A&A, 473, 

129 

Carretta, E., Bragaglia, A., Gratton, R. G., & Tosi, M. 2004, 

A&A, 422, 951 
— . 2005, A&A, 441, 131 

Casagrande, L., Schonrich, R., Asplund, M., Cassisi, S., Ramirez, 
I., Melendez, J., Bensby, T., & Feltzing, S. 2011, A&A, 530, 
A138-I- 

Cescutti, G., Matteucci, F., Frangois, P., & Chiappini, C. 2007, 
A&A, 462, 943 



& Kawata, D. 2005, 
., & Freeman, K. C. 2004, 
., Brook, C. B., & 
Pont, F. 2001, A&A, 372, 



Chabrier, G. 2001, ApJ, 554, 1274 

Chen, L., Hou, J. L., & Wang, J. J. 2003, AJ, 125, 1397 
Chiappini, C, Matteucci, F., & Gratton, R. 1997, ApJ, 477, 765 
Chiappini, C, Matteucci, F., & Romano, D. 2001, ApJ, 554, 1044 
Chiba, M., & Beers, T. C. 2000, AJ, 119, 2843 
Cresci, G., Mannucci, F., Maiolino, R., Marconi, A., Gnerucci, A., 

& Magrini, L. 2010, Nature, 467, 811 
Dalcanton, J. J., & Bernstein, R. A. 2002, AJ, 124, 1328 
de Jong, J. T. A., Yanny, B., Rix, H., Dolphin, A. E., Martin, 

N. F., & Beers, T. C. 2010, ApJ, 714, 663 
Dekel, A., Sari, R., & Ceverino, D. 2009, ApJ, 703, 785 
Di Matteo, P., Lehnert, M. D., Qu, Y., & van Driel, W. 2011, 

A&A, 525, L3-I- 

Dierickx, M., Klement, R., Rix, H.-W., & Liu, C. 2010, ApJ, 725, 
L186 

Dotter, A., Chaboyer, B., Jevremovic, D., Kostov, V., Baron, E., 

& Ferguson, J. W. 2008, ApJS, 178, 89 
Eisenstein, D. J., et al. 2011, AJ, 142, 72 
Elmegreen, B. G., & Elmegreen, D. M. 2005, ApJ, 627, 632 
— . 2006, ApJ, 650, 644 

FriedU, D., Benz, W., & Kennicutt, R. 1994, ApJ, 430, L105 
Friel, E. D., Jacobson, H. R., & Pilachowski, C. A. 2005, AJ, 129, 
2725 

— . 2010, AJ, 139, 1942 

Friel, E. D., Janes, K. A., Tavarez, M., Scott, J., Katsanis, R., 
Lotz, J., Hong, L., & Miller, N. 2002, AJ, 124, 2693 

Frinchaboy, P. M., Marino, A. F., Villanova, S., Carraro, G., 
Majewski, S. R., & Geisler, D. 2008, MNRAS, 391, 39 

Fukugita, M., Ichikawa, T., Gunn, J. E., Doi, M., Shimasaku, K., 
& Schneider, D. P. 1996, AJ, 111, 1748 



22 



Cheng et al. 



Gilmorc, G., & Rcid, N. 1983, MNRAS, 202, 1025 

Gilmorc, G., Wyse, R. F. G., & Jones, J. B. 1995, AJ, 109, 1095 

Girardi, L., Grebel, E. K., Odcnkirchcn, M., & Chiosi, C. 2004, 

A&A, 422, 205 
Goetz, M., & Koeppen, J. 1992, A&A, 262, 455 
Gunn, J. E., et al. 1998, AJ, 116, 3040 
— . 2006, A J, 131, 2332 

Hayashi, H., & Chiba, M. 2006, PASJ, 58, 835 
Haywood, M. 2008, MNRAS, 388, 1175 

Ibata, R., Irwin, M., Lewis, G., Ferguson, A. M. N., & Tanvir, N. 

2001, Nature, 412, 49 
Ivezic, Z., et al. 2008, ApJ, 684, 287 

Jacobson, H. R., Friel, E. D., & Pilachowski, C. A. 2008, AJ, 135, 
2341 

— . 2009, AJ, 137, 4753 

— . 2011, AJ, 141, 58 

Juric, M., et al. 2008, ApJ, 673, 864 

Kazantzidis, S., Bullock, J. S., Zcntncr, A. R., Kravtsov, A. V., & 

Moustakas, L. A. 2008, ApJ, 688, 254 
Kazantzidis, S., Zentner, A. R., Kravtsov, A. V., Bullock, J. S., & 

Debattista, V. R 2009, ApJ, 700, 1896 
Lacey, C. G., & Fall, S. M. 1985, ApJ, 290, 154 
Larson, R. B. 1976, MNRAS, 176, 31 
Lee, Y. S., ct al. 2008a, AJ, 136, 2022 
— . 2008b, AJ, 136, 2050 
— . 2011a, ApJ, 738, 187 
— . 2011b, AJ, 141, 90 

Lemasle, B., Frangois, P., Piersimoni, A., Pedicelli, S., Bono, G., 
Laney, 0. D., Primas, F., & Romaniello, M. 2008, A&A, 490, 
613 

Loebman, S. R., Roskar, R., Debattista, V. P., Ivezic, Z., Quinn, 

T. R., & Wadsley, J. 2011, ApJ, 737, 8 
Luck, R. E., Gicrcn, W. P., Andricvsky, S. M., Kovtyukh, V. V., 

Fouque, P., Pont, F., & Kienzle, F. 2003, A&A, 401, 939 
Luck, R. E., Kovtyukh, V. V., & Andrievsky, S. M. 2006, AJ, 

132, 902 

Luck, R. E., & Lambert, D. L. 2011, AJ, 142, 136 
Maciel, W. J., Lago, L. G., & Costa, R. D. D. 2005, A&A, 433, 
127 

Magrini, L., Randich, S., Zoccali, M., Jilkova, L., Carraro, G., 
Galli, D., Maiorca, E., & Busso, M. 2010, A&A, 523, A11+ 

Magrini, L., Sestito, P., Randich, S., & Galli, D. 2009, A&A, 494, 
95 

Martin, P., & Roy, J.-R. 1994, ApJ, 424, 599 
Martmez-Delgado, D., et al. 2010, AJ, 140, 962 
Matteucci, F., & Francois, P. 1989, MNRAS, 239, 885 
Minchev, I., & Famaey, B. 2010, ApJ, 722, 112 
Newberg, H. J., et al. 2002, ApJ, 569, 245 
Nordstrom, B., et al. 2004, A&A, 418, 989 
Pedicelli, S., et al. 2009, A&A, 504, 81 
Portinari, L., & Chiosi, C. 2000, A&A, 355, 929 
Prantzos, N., & Boissier, S. 2000, MNRAS, 313, 338 



Purcell, C. W., Kazantzidis, S., & Bullock, J. S. 2009, ApJ, 694, 
L98 

Randich, S., Sestito, P., Primas, F., Pallavicini, R., & Pasquini, 

L. 2006, A&A, 450, 557 
Read, J. I., Lake, G., Agertz, O., & Debattista, V. P. 2008, 

MNRAS, 389, 1041 
Reid, I. N., Gizis, J. E., & Hawley, S. L. 2002, AJ, 124, 2721 
Roskar, R., Debattista, V. P., Quinn, T. R., Stinson, G. S., & 

Wadsley, J. 2008a, ApJ, 684, L79 
Roskar, R., Debattista, V. P., Stinson, G. S., Quinn, T. R., 

Kaufmann, T., & Wadsley, J. 2008b, ApJ, 675, L65 
Sales, L. V., et al. 2009, MNRAS, 400, L61 

Sanchez-Blazquez, P., Courty, S., Gibson, B. K., & Brook, C. B. 

2009, MNRAS, 398, 591 
Schlegel, D. J., Finkbeiner, D. P., & Davis, M. 1998, ApJ, 500, 

525 

Schonrich, R., & Binney, J. 2009a, MNRAS, 396, 203 
— . 2009b, MNRAS, 399, 1145 

Sellwood, J. A., & Binney, J. J. 2002, MNRAS, 336, 785 
Sestito, P., Bragaglia, A., Randich, S., Carretta, E., Prisinzano, 

L., & Tosi, M. 2006, A&A, 458, 121 
Sestito, P., Bragaglia, A., Randich, S., Pallavicini, R., 

Andrievsky, S. M., & Korotin, S. A. 2008, A&A, 488, 943 
Sestito, P., Randich, S., & Bragaglia, A. 2007, A&A, 465, 185 
SmoUnski, J. P., et al. 2011, AJ, 141, 89 
Spitoni, E., & Matteucci, F. 2011, A&A, 531, A72+ 
Stanghellini, L., & Haywood, M. 2010, ApJ, 714, 1096 
Steinmetz, M., et al. 2006, AJ, 132, 1645 

Stewart, K. R., Bullock, J. S., Wechsler, R. H., Mailer, A. H., & 

Zentner, A. R. 2008, ApJ, 683, 597 
Stinson, G. S., Bailin, J., Couchman, H., Wadsley, J., Shen, S., 

Nickcrson, S., Brook, C, & Quinn, T. 2010, MNRAS, 408, 812 
Twarog, B. A., Ashman, K. M., & Anthony-Twarog, B. J. 1997, 

AJ, 114, 2556 

Venn, K. A., Irwin, M., Shetrone, M. D., Tout, C. A., Hill, V., & 

Tolstoy, E. 2004, AJ, 128, 1177 
Villalobos, A., & Helmi, A. 2008, MNRAS, 391, 1806 
Villanova, S., Carraro, G., Bresolin, F., & Patat, F. 2005, AJ, 

130, 652 

Wilson, M. L., ct al. 2011, MNRAS, 413, 2235 

Yanny, B., et al. 2009, AJ, 137, 4377 

Yoachim, P., & Dalcanton, J. J. 2005, ApJ, 624, 701 

— . 2006, AJ, 131, 226 

— . 2008, ApJ, 683, 707 

Yong, D., Carney, B. W., & Teixera de Almeida, M. L. 2005, AJ, 

130, 597 

Yong, D., Carney, B. W., Teixera de Almeida, M. L., & Pohl, 

B. L. 2006, AJ, 131, 2256 
York, D. G., et al. 2000, AJ, 120, 1579 
Yoshii, Y. 1982, PASJ, 34, 365 



